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A PREDICTIVE  SCHEME  FOR  THE  BLAST  ENVIRONMENT  OF  ARMY  WEAPONS 


PREFACE 


An  Important  Note  to  the  Reader 

This  memorandum  report  is  the  first  of  three  parts  of  an  analysis 
which  predicts  the  blast  wave  properties  produced  by  US  Army  weapons. 
Experience  and  reflection  have  ultimately  shown  a better  way  to  perform 
the  analysis  in  this  as  in  many  other  endeavors:  Thus,  a three  part 
report . 

The  preliminary  analysis  in  this  report  was  designed  to  establish 
the  validity  of  a computationally  simple  approach  to  the  problem. 

Emphasis  was  placed  on  elucidation  of  the  theories  with  minimal  compari- 
son to  data.  Further  work,  to  be  presented  in  Part  II,  has  applied  the 
theory  to  a wide  range  of  calibers.  These  extensions  '^clean  up”  some 
areas  which  have  not  previously  been  founded  in  rigorous  analyses. 

This  new  theory  has  not  been  completely  exploited;  e.g.,  work  is  in 
progress  to  apply  the  theory  for  prediction  of  the  rearward  blast  field 
produced  by  recoilless  rifles.  This  work  will  appear  as  Part  III. 

The  analysis  was  initiated  as  a result  of  the  promulgation  of  MIL- 
STD-1474A(MI) , "Noise  Limits  for  Army  Materiel".  Where  previous  standards 
had  provided  no  serious  constraints  on  the  number  of  training  rounds 
that  could  be  fired  by  crews  serving  artillery  pieces,  the  new  standard 
introduced  requirements  which  might  make  major  changes  in  training  sched- 
ules --  and  costs.  In  particular,  the  section  on  impulse  noise  for 
personnel  - occupied  areas  limit  "blast^^  exposures  ranging  from  IQOO  to 
no  exposures  per  day. 

The  potential  cost  of  the  experimental  program  to  certify  that  all 
service  locations  for  crew  served  weapons  became  a matter  of  concern  to 
the  Director  of  AMSAA.  In  turn,  he  suggested  that  experimental  work 
could  be  minimized  if  the  BRL  produced  an  a phlohl  predictive  theory  to 
guide  such  blast  level  experiments  as  would  be  necessary  to  implement  the 
new  Mil  - Standard.  Specifically,  it  was  suggested  that  BRL  Report  1019 
had  experimental  data  that  could  provide  guidance  in  formulating  a theory. 
In  addition,  the  smaller  computational  scheme  was  to  be  preferred  over 
hydro-codes  since  only  scalar  features  of  the  flow  field  are  involved  in 
determination  of  compliance  with  the  Mil-Standard. 

To  the  individual  who  feels  our  approach  is  a fortuitous  process, 
we  offer  the  defense  of  H,  Bethe:  "If  this  is  not  the  correct  theory,  it 
is  still  an  excellent  way  to  correlate  the  data". 
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This  report  presents  our  application  of  existing  technology  to 
the  problem.  We  are  indebted  to  our  colleagues  in  BRL,  AMSAA,  and  HEL 
for  their  cooperation,  both  in  providing  data  and  for  the  counsel  and 
advice  we  received.  Foremost,  however,  we  are  indebted  to  Francis 
Porzel,  without  whose  theories  this  computational  scheme  would  not  have 
been  possible. 
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INTRODUCTION 


This  report  presents  a method  for  predicting  the  muzzle  blast 
overpressures  and  pressure  pulse  length  emanating  from  guns  during 
firing.  The  study  was  undertaken  to  determine  if  ear  protection  re- 
quired for  various  guns,  can  be  predicted  a priori ^ i.e.,  without  the 
necessity  of  heretofore  required  extensive  experimentation.  Our  objec- 
tive is  to  define  a method  which  is  independent  of  overpressure  and 
pulse  length  data  (i.e.,  no  curve  fitting)  and  which  is  simple  and 
economical  to  use.  The  initial  investigation  is  limited  to  large  cali- 
ber guns  without  muzzle  devices. 

Previous  investigations  fall  into  two  categories:  (1)  Large 
computer  programs  (the  so-called  hydrocodes)  and  (2)  massive  correla- 
tions of  data.  Representative  examples  of  the  first  category  include 
the  works  of  Schmidt^ Zoltani^  and  Ranlet^.  These  investigations 
provide  insight  into  the  basic  physical  processes  but  they  require 
extensive  computations  to  produce  a solution.  An  excellent  investiga- 
tion within  the  second  category  is  the  work  of  Westline^  which  estimates 
the  blast  overpressure  for  various  guns.  His  results  include  empirical 
constants  which  are  necessary  in  describing  the  energy  available  to  the 
blast  from  a specific  gun.  Neither  the  computational  nor  the  empirical 
approach  is  acceptable  in  view  of  the  obj ective- -simple,  general  and. 
economical  predictive  capability. 

The  most  extensive  work  in  blast  overpressure  prediction  has  been 
in  the  area  of  blast  from  explosions.^  By  this  we  mean  that  area  in 
which  spherical  symmetry  is  obeyed.  The  technical  literature  of  blast 
(explosions)  abound  with  similarity  parameters,  scaling  laws,  governing 


Schmidt^  R.E.  Shear ^ ^^The  Flow  Field  About  the  Muzzle  of  an  M16 
Rifle,  BRL  Report  No.  1692,  Jan.,  1974  (AD  #9166461) 

^E.M.  Schmidt,  R.E.  Shear,  ^^Launoh  Dynamics  of  a Single  Flachette  Round, 
BRL  Report  No.  1810,  Aug.,  1975.  (AD  #B006781L) 

^C.K.  Zoultani,  ^^Evaluation  of  the  Computer  Codes  BLAST,  DORF,  HELP  and 
HEMP  for  Suitability  of  Under  expanded  Jet  Flow  Calculations,  " BRL  Report 
No.  1659,  Aug.,  1973  (AD768708). 

^J.  Ranlet,  J . Erdos,  ^^Muzzle  Blast  Flow  Field  Calculations, ''  BRL  Con-^ 
tract  Report  No.  297,  Apr.,  1976.  (AD  #B011967L) 

^P.S.  Westine,  J.C.  Hokanson,  ^^Prediction  of  Stand-off  Distances  to 
Prevent  Loss  of  Hearing  from  Muzzle  Blast, " Rock  Island  Arsenal  Report 
No.  R-CR-75--003,  Feb.,  1975  (AD/A-005274) . 

^N.E.  Baker,  ^^Explosions  in  Air,^^  University  of  Texas  Press,  Austin  and 
London,  1973. 
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relations  for  weak  and  strong  shock  waves,  etc.  We  have  applied  this 
extensive  work  to  the  gun  blast-field  problem.  This  was  done  because 
at  locations  far  from  the  source  of  the  explosion,  i.e.,  in  the  weak 
shock  regime,  asymmetrical  effects  are  washed  out  and  the  shock  wave 
can  be  modeled  as  though  it  originated  from  a spherical  source  of 
finite  radius. 

The  gun  blast  problem  requires  analyses  in  addition  to  spherical 
explosion  theory.  The  first  is  an  estimate  for  the  equivalent  explo- 
sive yield  which  is  produced  by  the  gases  escaping  through  the  muzzle. 
Using  the  propellant  charge  energy  (reduced  by  the  projectile  energy) 
is  not  correct  because  there  are  losses  in  the  interior  of  the  gun 
(heating,  recoil,  boundary  layer  generation,  etc.).  The  second  analysis 
addresses  the  (highly  directional)  energy  release  from  the  gun;  which 
is  now  more  cylindrically  than  spherically  shaped.  Lastly,  other 
theories  do  not  admit  determination  of  pulse  length  without  recourse  to 
the  above  mentioned  correlation  techniques.  The  principal  contribution 
of  this  report  is  a general  treatment  of  these  areas. 

Assuming  that  the  above  problems  have  been  solved,  we  still  require 
a theory  for  the  spherical  blast  wave.  Scaling  laws  provide  no  help  in 
an  a priori  analysis.  A similitude  analysis  is  required  and  perhaps  the 
most  well  known  is  that  by  Sir  G.I.  Taylor  in  1950^.  His  analysis 
reduced  the  several  governing  differential  equations  to  one  which  can  be 
solved  by  numerical  quadratures.  Seeking  an  even  simpler  form  (mathe- 
matically) we  have  selected  the  Unified  Theory  of  Explosions  (UTE)  as 
advanced  by  F.B.  Porzel®,  UTE  is  a comprehensive  theory  for  explosions 
which  offers  simple  analytic  expressions  for  predicting  blast  parameters. 
This  theory  has  been  applied  to  a variety  of  explosions,  from  thermonu- 
clear to  small  HE  charges,  and  found  to  be  valid  to  within  a few  percent. 

The  following  sections  include:  II.  Interior  losses;  III.  Theory; 
IV.  Presentation  of  Data,  Comparison  between  theory  and  experiment, 
areas  not  addressed  by  our  analysis,  and  an  Appendix  describing  a 
computer  code  listing  and  sample  output. 


Taylor y ^^The  Formation  of  a Blast  Wave  by  a Very  Intense  Explosion: 
I Theoretical  Discussion^  Froc.  R.  Soc.  A.  201^  259-174^  (1950). 


Q 

F.B.  Porzel,  "Introduction  to  a Unified  Theory  of  Explosions  (UTE)," 
NOLTR-7 2-209,  Sept.,  1972  (AD-758000). 
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II.  THEORY 


Analysis  of  the  muzzle  blast  must  be  preceded  by  analyses  (or  at 
least  definition)  of  the  processes  that  occur  in  the  gun  interior,  at 
the  muzzle  and  of  those  which  contribute  to  the  shock  wave  propagation. 
The  interior  processes  include  compression  and  expansion  waves,  bound- 
ary layer  flows,  complex  chemical  reactions,  heat  transfer  from  the 
fluid  to  solid,  friction,  etc.  The  muzzle  flow  field  develops  with 
many  similarities  to  the  processes  of  underexpanded  nozzles.  There  is 
a complicated  interaction  between  the  fluid  of  the  gun  and  the  atmos- 
phere with  the  formation  of  a so-called  shock  bottle  and  associated 
Mach  disks.  Eventually,  at  a distance  from  the  muzzle  the  flow  field 
becomes  more  regular,  resembling  a geometrically  transformed  version 
of  the  blast  wave  that  would  be  generated  by  .an  explosion.  Current 
physical  understanding  of  spatial  non-uniformity  is  that  the  greatest 
transport  of  any  quantity  is  in  the  direction  of  the  largest  gradient; 
hence,  as  the  asymmetrical  wave  propagates  outward  it  should  approach 
spherical  symmetry. 

We  propose  that  the  complicated  muzzle  blast  phenomena  may  be 
modeled  as  if  it  were  initiated  by  a spherical  explosive  charge. . 

Ear  protection  required  by  military  regulations  is  given  in  terms 
of  the  peak  pressure  and  duration  of  the  shock  wave;  that  is,  in  terms 
of  the  energy  in  the  shock  wave.  Consequently,  our  analysis  is  an 
account  of  the  way  in  which  the  initial  propellant  energy  is  partitioned; 
into  heat,  friction,  internal  energy,  propellant  and  projectile  kinetic 
energy,  etc.  Finally,  it  is  an  account  of  that  energy  available  to  the 
shock  wave.  As  we  shall  show,  it  is  sufficient  in  the  muzzle  blast 
problem  to  monitor  the  energy  transport  processes--we  do  not  require 
a detailed  knowledge  of  the  flow. 

The  assumption  of  an  equivalency  between  the  muzzle  blast  and 
spherical  explosions  has  an  important  consequence.  With  this  assump- 
tion we  are  able  to  use  all  of  the  analyses  which  occurred  during  the 
early  phases  of  the  nuclear  explosion  era  for  the  solution  of  the  muz- 
zle blast  problem.  An  excellent  exposition  of  many  of  the  theories  and 
experimental  data  can  be  found  in  Baker^ . Curiously  missing  are  the 
theories  we  shall  use--those  of  F.B.  Porzel . In  the  following  sections 
we  shall  derive  and  demonstrate  Porzel *s  theories  and  show  how  these 
analyses  can  be  applied  to  the  muzzle  blast  overpressure  problem.* 

An  important,  underlying  notion  must  be  addressed  first  - the 
division  of  energies  for  shock  wave  overpressure  prediction.  In  the 
following  sections  we  shall  speak  of  prompt  energy  and  waste  heat. 

Prompt  energy  is  that  energy  which  is  (promptly)  available  for  driving 

*Most  of  PorzeT^s  work  is  available  in  institutional  or  corporate 
documents.  But;  since  Porzel *s  work  is  neither  available  in  archival 
journals  nor  included  in  surveys  like  that  of  Baker,  the  authors  have 
chosen  to  include  their  derivation  of  Porzel *s  work  here. 
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the  shock  wave.  Any  energy  which  is  lost  from  the  shock  wave,  by 
whatever  means,  is  termed  waste  heat.  We  emphasize  that  the  waste 
heat  is  not  necessarily  lost  from  the  system;  it  simply  is  not  avail- 
able to  the  shock  wave.  This  differentiation  is  important  because  in 
many  instances  one  would  be  led  to  say  that  neglecting  this  energy  or  ; 
that  energy  will  have  severe  consequences  when  predicting  certain 
phenomena  other  than  the  shock  wave.  We  are  concerned  with  only  the 
shock  wave,  its  strength  and  duration.  Any  energy  which  does  not  i 
support  the  wave  is,  for  our  purposes,  wasted  (or  delayed)*. 

Prompt  energy  includes  the  pressure  volume  work  in  an  expansion 
process  and  the  kinetic  energy  of  ordered  motion  imparted  to  the  gas 
during  the  expansion.  All  other  energy  is  waste  energy.  We  can  study 
the  waste  heat  graphically  with  the  help  of  the  P-v  diagram  for  an 
adiabatic  expansion  process  (Figure  1) . 


FIGURE  1 

From  the  Rankine-Hugoniot  relations^  , we  find  that  the  rectangle 
bounded  by  P = constant,  v = constant  and  P = 0,  Va  = constant  is 
the  total  energy: 


ej  = P - v). 

The  upper  right  triangle  in  the  rectangle  is  the  kinetic  energy: 


The  remaining  trapazoidal  area  is  the  internal  energy: 

i O’  * 


*Delayed  energy  is  a subset  of  waste  heat. 

^Y.B.  Zel ' doTviah,  Yu  P.  Raizev,  "Physics  of  Shock  Waves  and  High- 
Temperature  Hydrodynamic  Phenomena^  " Academic  Press^  Hew  York  and 
London,  1966 
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In  an  adiabatic  expansion  process  (where  real  gas  effects,  such  as 
ionization  etc.  are  neglected)  the  gas  pressure-volume  will  generally 
follow  the  curved  line,  the  so-called  Hugoniot  or  shock  adiabat.  The 
area  between  the  adiabat  and  the  straight  line  connecting  Pv  and  PaVa 
is  the  waste  heat.  At  an  infinite  distance  from  the  source  of  the 
blast  all  of  the  prompt  energy  has  been  converted  to  waste  heat. 

We  can  estimate  the  ratio  of  production  of  waste  heat  to  the 
total  energy  using  the  Rankine-Hugoniot  relations  for  changes  across 
a shock  wave : 


P+P 


a 


2P 


For  P » Pa  the  internal  energy  is  approximately  1/2  the  total,  hence, 
half  the  energy  is  subject  to  waste  (the  kinetic  energy  is  not  subject 
to  waste)  . In  the  acoustic  approximation,  P 'v  P^  and  all  the  energy 
is  wasted. 

The  following  solution  scheme  is  separated  into  two  distinct 
analyses.  The  first  is  concerned  entirely  with  the  estimate  for 
the  prompt  energy  as  the  gases  leave  the  muzzle.  These  calculations 
are  used  to  predict  a yield  and  radius  for  an  equivalent  spherical 
explosion  which  is  treated  as  a separate  problem. 
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A.  INTERIOR  LOSSES 


Losses  interior  to  the  gun  (keeping  in  mind  the  objective  of  finding 
values  for  waste  heat  and  prompt  energy)  have  been  analyzed  as: 

(1)  Kinetic  energy  given  to  the  projectile  and  so  lost  to  the 
supply  of  prompt  energy;  (2)  other  interior  losses  or  delays  of  energy 
which  amount  to  about  5/6 's  of  the  otherwise  available  energy*; 

(3)  energy  loss  due  to  expansion;  and,  (4)  energy  loss  due  to  tube 
roughness  (in  this  case,  rifling). 

Item  (1)  is  self  evident.  Item  (2)  is  compatible  with  the  mass 
effect  (Section  B-2)  which  states  that  initially  the  shock  wave  is 
driven  by  the  mass  of  the  explosive.  If  we  view  the  propellant  as  a 
cube  (Figure  2)  we  argue  that  only  1/6  of  the  mass  is  directed  so  as 
to  be  promptly  available.  The  remaining  mass  which  travels  laterally 
and  to  the  rear  eventually  exits  the  gun  but  because  it  is  delayed, 
it  does  not  directly  contribute  to  the  exiting  shock  wave. 


RECOIL 

DIRECTION 


BARREL  AXIS 


Regarding  Item  (3) , we  defined  the  prompt  energy  as  including  the 
kinetic  energy  of  ordered  motion.  When  a projectile  travels  down  a 
barrel  we  expect  a turbulent  boundary  layer  (in  which  the  ordered  motion 
becomes  random)  to  be  formed,.  After  a sufficient  length  of  the  barrel 
has  been  traversed  we  further  expect  the  boundary  layer  to  close  on 
itself;  that  is,  at  some  distance  behind  the  projectile  there  will  be 
a point  at  which  the  turbulent  motion  extends  across  the  entire  diameter 
of  the  barrel.  We  treat  the  energy  in  the  random  motion  as  delayed 
and  hence,  wasted. 


*This  is  consistent  with  Porzel ' s work  applied  to  conical  shock  tubes, 
lO 

, F.B.  Porzet^  ^^Correlation  of  Blast  Simulators  with  a Unified  Theory 
of  Explosions ^ ■*  3rd  International  Symposium  on  Military  Application  of 
Blast  Simulators^  Schwetzinger^  Germany^  Sept. j 1972. 
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Elementary  flow  theory  tells  us  that  when  the  velocity  of  a gas 
increases  the  static  pressure  decreases  according  to 


.2 


P « K - aV 


That  is,  given  the  same  reservoir  conditions  (unchanged  K) , higher 
velocity  gases  have  lower  static  pressures.  The  strength  of  the  shock 
wave  which  is  formed  by  the  escaping  gases  will  be  determined  by  the 
static  pressure  at  the  projectile  base  (relative  to  an  inertial  refer- 
ence frame) . If  we  apply  these  notions  to  the  material  velocity 
(velocity  of  ordered  motion)  of  the  flow  we  can  view  the  flow  as  pro- 
ceeding from  some  stagnation  pressure  at  the  point  of  closure 
(’^material  0)  to  some  static  pressure  at  the  projectile  base.  This 
represents  an  energy  reduction  (via  the  pressure  expansion)  and  is 
illustrated  in  Figure  3. 


/////////////////////// 


7//////////////////  / / // 


FIGURE  3 


The  last  significant  loss  is  due  to  the  roughness  (rifling)  in  the 
barrel.  The  energy  in  the  fluid  which  interacts  with  (and  is  trapped 
by)  the  rough  wall  is  delayed  and  is  not  available. 

Item  (1)  is  easily  determined  from  the  muzzle  velocity  and  mass  of 
the  projectile.  Item  (2),  which  is  not  universally  justifiable,  is 
reasonable  within  the  context  of  our  energy  definitions.  The  "proof 
of  the  puddin"  is  in  the  comparison  with  experiment  (Section  III). 
Following  is  the  approach  for  determining  the  losses  caused  by  effects 
(3)  and  (4).  These  analyses  are  obtained  from  Porzel^^. 

A. 1 Energy  Loss  in  the  Presence  of  Boundary  Layer  Choke 

We  are  concerned  with  a loss  of  (stagnation)  energy  in  the  presence 
of  a turbulent  boundary  layer  in  the  gun  barrel.  The  boundary  layer, 
created  by  the  moving  projectile,  grows  as  the  projectile  travels  the 
length  of  the  tube,  but  at  some  point  the  boundary  layers  created  on 
the  sides  meet  at  the  axis  of  the  barrel.  We  shall  speak  of  a closure 
of  the  boundary  layer  (Figure  3) . Once  this  closure  has  occurred  the 
energy  of  the  gas  is  definitely  separated  into  a region  of  ordered  motion 
and  a region  where  the  significant  portion  of  the  energy  is  in  the  ran- 
dom motion  of  the  fluid  particles. 


11 

F.B,  Porzelj  "Study  of  Shook  Irr^edance  Ef feats  in  a Rough  Waited 
Tunnel,  " Institute  for  Defense  Analysis  Researah  Paper  P-330, 

Mar.,  1969  (AD684790). 
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Once  closure  has  occurred  we  analyze  the  flow  as  a Bernoulli 
expansion  from  the  stagnation  pressure  at  the  closure  (the  material 
velocity  or  ordered  motion  is  essentially  zero)  to  the  static  pressure 
at  the  base  of  the  projectile.  We  associate  the  stagnation  pressure 
with  the  maximum  pressure  and  associate  the  pressure  at  the  projectile 
base  with  the  initial  overpressure  of  the  shock  wave.  To  make  the 
analysis  tractable  we  assume  there  are  no  compression  or  expansion 
waves  in  the  region  from  the  closure  point  to  the  projectile  (i.e.,  we 
can  treat  the  expansion  as  adiabatic) . 

Our  only  interest  is  fluid  motion  along  the  axis  of  the  barrel; 
hence,  we  can  write  the  equation  describing  conservation  of  momentum 
in  a moving  coordinate  system  attached  to  the  projectile  as 


9^ 

3z 


3u 


— = 0 
3t 


(1) 


in  which  the  material  velocity  is  represented  by  u. 

One  may  assume  that  once  the  boundary  layer  closure  has  occurred 
the  distance  from  the  closure  point  to  the  base  of  the  projectile 
remains  constant.  That  is,  steady  state  prevails  and  3u/3t  = 0.  (In 
Section  IIl.A.l  we  show  that  the  ordered  velocity  is  of  the  same  order 
as  the  random  speed  in  one  direction.)  Additionally,  since  the  flow 
is  adiabatic. 


P = Ap^ 


(2) 


Substituting  Equation  (2)  into  (1) , we  obtain 


which  can  be  integrated  to  give 


pi-i/r  , 1 ,2 

1--  ^ 

Y 


= constant 
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Since 


1/y 


P 


and,  the  speed  of  sound,  a,  is  given  by 


1/2 

a = 


we  finally  obtain 


constant. 


(3) 


Equation  (3)  relates  the  material  velocity  and  speed  of  sound  at 
the  closure  point  to  the  same  quantities  at  the  base  of  the  projectile. 
This  equation  applies  to  the  barrel  reference  frame  provided  the  assump- 
tions concerning  constant  closure-to-proj ectile-base  distance  and  isen- 
tropic  flow  are  not  violated. 

The  constant  in  Equation  (3)  can  be  evaluated  at  the  closure  point 
where  u « 1,  u "v  0.  Thus  Equation  (3)  becomes 


(4) 


where  the  subscript  c denotes  the  choke  point. 


From  Equation  (4)  we  obtain 


(5) 


The  adiabatic  relation  allows  us  to  write 
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which  can  be  reduced  to 


Combining  Equations  (5)  and  (6)  produce  Equation  (7) : 


Y-1 


1 


+ 


111 

2 


(6) 


(7) 


Let  us  assume,  for  the  moment,  that  the  length  of  the  barrel  is 
equal  to  the  distance  from  the  breech  to  the  projectile  base  at  the 
moment  of  boundary  layer  closure.  The  projectile  would  rapidly  leave 
the  driving  gases  because  of  the  lateral  expansion  of  these  gases. 

The  ordered  motion  is  supersonic  with  respect  to  the  ambient  air 
hence  a shock  wave  will  form  at  the  gas  leading  edge.  At  that  instant 
we  can  use  the  Rankine-Hugoniot  relations  to  determine  the  pressure 
ratio  across  this  shock.  (The  shock  wave  which  in  fact  produces  the 
strongest  overpressure) . Specifically,  we  find^ 


and  Equation  (7)  becomes 

Y-1 


or,  denoting  the  pressure  energy  driving  the  shock  wave  by  Prs^ 


Equation  (8)  relates  the  reduction  in  the  energy  from  the  (essen- 
tially) initial  pressure  ratio,  , to  the  pressure  ratio  energy 
promptly  available  to  the  shock  wave,  P;j^  . 
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For  current  guns  of  interest,  pressure  ratios  in  excess  of  100  are 
not  unusual.  Neglecting  numbers  of  order  unity  (compared  with  100)  and 
using  a specific  heat  ratio  of  1.25  Equation  (8)  becomes,  approximately. 


This  means  that  the  effect  of  the  turbulent  choke  is  a reduction  in  the 
available  energy  for  driving  the  shock  wave  by  a factor  of  approximately 
19. 

In  order  to  estimate  the  point  at  which  closure  occurs  we  note  that 
experimental  results  involving  shock  tubes  with  walls  of  known  roughness 
give  (denoting  a roughness  factor  by  H) : 


L 15 
D-rO.! 


(9) 


as  a very  reasonable  fit  of  the  data^^.  The  roughness  factor  is  the 
ratio  of  the  roughness  height,  h,  to  the  unimpeded  diameter  of  the  tube. 

A. 2 Energy  Loss  Due  to  Impedance  to  the  Flow  by  the  Rough  Wall 

The  prompt  energy  includes  the  ordered  kinetic  energy  of  the  flow. 
Near  the  wall  where  the  roughness  of  the  wall  can  be  felt  by  the  flow, 
there  is  a local  decrease  in  the  kinetic  energy  as  the  flow  encounters 
protuberances.  These  energy  losses  occur  in  addition  to  the  boundary 
layer  and  we  expect  them  to  appear  near  the  base  of  the  projectile. 

We  shall  assume  that  the  change  in  the  total  energy  is  proportional 
to  the  kinetic  energy  (per  unit  volume),  and  the  volume  subtended 

by  the  average  roughness  of  the  barrel,  which  irselt  is  a product  of 
the  roughness  height,  h,  perimeter  S and  distance  dL.  That  is 


dE^  = -a  Cj^gShdL 


(10) 


where  the  constant  of  proportionality,  a,  can  be  viewed  as  an  absorp- 
tion coefficient. 

It  is  reasonable  to  assume  that  the  losses  occur  in  a volume  of 
dimensions  A*D,  where  A is  the  cross-sectional  area  and  D is  the 
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diameter.  If  the  total  energy  per  unit  volume  is  denoted  by  ej,  then 
Equation  (10)  can  be  written 


or 


A-Dde^  = - a e^SHdL 
1 Kc 


de^  = - 


a e 


S h 
KE  A D 


dL 


(11) 


We  can  compare  this  relation  with  the  classical  exponential  decay 
law  by  writing  it  in  the  form 


S h 

“ad 


e 


T ■ 


That  is, 


de 


T 


dL 


CC 


(constant)  • fn(P) 


e 


T 


where  the  absorption  coefficient  is  a pressure  dependent  function. 

We  can  write  Equation  (11)  in  dimensionless  form  by  setting  the 
diameter,  D,  equal  to  the  hydraulic  diameter,  i.e.. 


S _ 4 
A - D ■ 


Equation  (11)  becomes 


If  we  define 


L 

d 
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X 


we  obtain 


de„ 


4 ot  H dx  . 


TE 


From  the  Rankine-Hugoniot  relations,  it  is  known  that 


and 


(12) 


V yP  +1 
o r 


P +y 
r 


• '■r  ' '’/'’a-  " ' 


hence , 


Et  - 


P +y 
r 


The  ratio  of  change  of  eT  with  the  pressure  ratio,  P^,  is 


(y-l)P  (P  +2yP  -y) 
de  = V ^ dP 


03) 


Again,  beginning  with  the  Rankine-Hugoniot  relations  we  can  find  that 

,2 


'KE 


(P-l)Pa(V^^‘ 

2(P^+y) 


(14) 


thus,  the  left  hand  side  cfuSquation  (12)  is 


d^T  +2yP  -y) 

— - = 2 — 2 dP 

''kE  (P+y)(P-l)^ 


(15) 
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If  we  define 


P-P 

AP  = -5-^ 
r P 


-=  P -1 
r 


and 


Equation  (15)  becomes 

de 


(AP  ) + 2e(AP  ) + 3 

— -=  2 ^ d(AP  ) . 

^KE  (AP  ) [(AP  ) + 3]  "" 


(16) 


Combining  Equations  (16)  and  (12)  the  incremental  energy  loss  due  to 
the  flow  impedance  is  described  by 


(AP^)  +23(AP^)+3 
(AP^)"(AP^+3) 


d(AP  ) = - 2 a 11  dx  . 


Integrating,  we  obtain 


— - — In  (AP  ) - In  (AP  + 3)  - 77:—  = constant  - 2 a H dx 
3 r 3 ^ r AP 

r 


Denoting  the  left  side  by  I,  this  relation  determines  the  energy  loss 
between  two  positions  (1  § 2)  down  the  barrel  in  the  form 


(17) 
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B.  THE  UNIFIED  THEORY  OF  EXPLOSIONS  - UTE 


Next,  we  shall  determine  the  shock  wave  properties  by  use  of  the 
Unified  Theory  of  Explosions  (UTE)  as  advanced  by  F.B.  Porzel®.  UTE 
is  a comprehensive  theory  providing  simple  analytic  expressions  for 
blast  parameter  determination.  UTE  in  total  covers  a gamut  of  condi- 
tions, geometries,  etc.;  however,  we  will  examine  only  those  parts 
pertinent  to  our  problem.  The  concepts  of  prompt  energy  and  waste 
heat  were  introduced  in  the  beginning  of  this  chapter.  We  now  wish 
to  show  how  this  division  of  energies  permits  calculation  of  the 
overpressure  for  various  ranges  from  the  muzzle. 

The  prompt  energy  was  defined  as  the  kinetic  energy  of  ordered 
motion  plus  the  static  overpressure  energy.  All  other  energy  which 
does  not  directly  support  the  wave  was  defined  as  waste  heat. 

Hence,  in  general  we  can  write 


(18) 


where  W is  the  pressure  volume  energy,  W = JPdv,  ®kE  kinetic 

energy  per  unit  mass  (dynamic  pressure)  and  Q is  the  waste  neat. 

The  total  prompt  energy  Y (R) , the  integral  of  the  prompt  energy, 
defined  by 


(19) 


o 


is  of  particular  interest.  The  boundary  conditions  for  the  shock 
expansion  process  are: 

(1)  At  R = R , the  initial  charge  radius,  Y,  is  the  hydrodynamic 
yield  of  the  explosion,  Y^;  and 

(2)  as  R approaches  infinity  the  shock  wave  must  be  completely 
dissipated  hence  Y (“)  ->  0. 

Substitution  of  Equation  (18)  into  (19)  gives 


o o 
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But  boundary  condition  (2)  requires  that 


j-TT  = 47r  / Qr^dr 


o 


hence 


Y(R)  = 47r  / Qr^dr 
R 


(20a) 


from  which  the  rate  o£  loss  o£  Y is  £ound  to  be 


(20b) 


Speci£i cation  o£  Q permits  determination  o£  the  prompt  energy  and  hence, 
o£  the  static  overpressure.  The  UTE  becomes  tractable  i£  Q can  be 
speci£ied. 

The  abstraction  which  makes  UTE  tractable  is  the  QZQ  hypothesis 
which  states; 


QZ^  = constant 


(21) 


where  Z is  a mass  corrected  radius.  This  relation  permits  Equation 
(20a)  to  be  integrated  in  closed  £orm  yielding  a simple  analytic  equa- 
tion £or  the  behavior  o£  the  prompt  energy.  Justi£ication  o£  this 
hypothesis  requires  development  o£  an  equation  o£  state,  determination 
o£  the  mass  corrected  radius  and  the  introductiofi  o£  a £orm  £actor. 

The  £ollowing  sections  develop  the  three  concepts  with  the  £inal  section 
devoted  to  proo£  o£  the  above  relation  (21) . 

B. 1 The  Generalized  Equation  o£  State  - GES 

Our  analysis  is  designed  to  cover  the  spectrum  o£  gaseous  states 
£rom  ambient  to  highly  compressed.  To  do  this,  we  require  an  equation 
o£  state  capable  o£  being  extended  into  the  dense  gaseous  state.  The 
classical  equation  o£  state  does  not  include  interactions  between  par- 
ticles caused  by  the  long  range  £orce.  It  is  derived  by  assuming  spheri- 
cal particles  without  an  interaction  potential. 
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Landau  and  Lifshitz^^  develop  a correction  term  to  the  perfect 
gas  law  by  assuming  the  gas  is  sufficiently  dense  that  binary  collisions 
are  important  but  triple,  quartic,  etc.,  collisions  may  be  neglected. 
Noting  that  the  pressure  can  be  found  from  the  Gibbs  free  energy,  F, 
by 


P = 


iE 

3v 


they  find  that,  for  binary  collisions,  the  free  energy  is  given  by 


F = Fp  + K(T)/v 


where  N is  the  total  number  of  particles,  Fp  is  the  free  energy  for  the 
perfect  gas  state,  v is  the  volume  and  K(T)  is  proportional  to  the  two 
particle  interaction  potential,  Ui2>  via 


K(T)  cc  / 


/kT 


dv 


We  can  generalize*  this  result  by  writing 


P = E 5 ./v’^i 
1 


where  the  and  are  constants  embodying  the  interparticle  inter- 
actions; these  constants  are  determined  from  the  thermodynamics  of  the 
processes.  We  can  write  this  relation  relative  to  the  ambient  pressure. 
Pa,  as 


a 


a 


which  may  be  reduced  to 


(22) 


*Alternatively  one  can  view  the  prior  results  as  a specialization  of 
Equation  22  (which  follows) . 

LandaUy  EM.  Lifshitz^  ^Statistical  Fhysios^^^  Pergamon  Press 
Ltd. y London j 1938 
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This  is  the  Generalized  Equation  of  State  (GES)  in  the  UTE  used  to 
calculate  the  waste  heat. 

Restricting  our  analyses  to  those  regimes  where  real  gas  effects 
can  be  neglected,  we  have 


P - 


One  constant  can  be  determined  by  recalling  that  a^^  = (dp/dp)  , i.e., 

Si 


These  equations  are  combined  to  yield 


Use  is  made  o£  two  characteristics  o£  a per£ect  gas  (with  constant 
ratio  o£  speci£ic  heats),  the  adiabatic  relation 


P = (constant)  x 


(23) 


and 


P 

a 


P a 
a a 


to  obtain 


2 

p a 

P - P = _a_a_ 
a Y 


(24) 


and  thus  we  conclude  that  on  the  average,  n = y- 
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B.2  The  Mass  Effect  (MEZ)  - The  Mass  Corrected  Radius  - Z 


Blast  energy  is  initially  contained  in  the  energetic  material  which 
produces  the  explosion,  i.e.,  the  propellant.  Since  no  propellant 
burns  instantaneously  or  (rarely)  to  completion,  at  initiation  the 
shock  wave  is  driven  by  the  products  of  the  reaction  and  by  the  as  yet 
unburned  particulate  matter.  As  the  shock  wave  expands  and  air  is 
engulfed,  some  of  the  energy  is  transferred  to  the  air  and  continues 
to  drive  the  shock.  Eventually,  the  shock  wave  leaves  the  residual 
mass  because  viscous  drag  on  the  particles  reduces  their  speed  rela- 
tive to  the  shock  speed.  The  effect  is  a reduction  in  the  pressure 
since,  as  in  the  case  of  smoke  which  has  a greater  specific  heat,  the 
energy  density  in  the  residual  mass  is  greater  than  that  which  would 
exist  if  only  air  were  present. 

In  MEZ,  the  energy  is  assumed  to  be  distributed  between  the  parti- 
culate mass  and  air  in  direct  proportion  to  their  relative  masses,  that 
is. 


BM 

prompt  energy  = (prompt  energy)  ^^^[1  + ^ ^ ] (25) 

3 


where  M is  the  particulate^mass , B is  the  ratio  of  the  mass  prompt 
energy  to  that  of  air  and  ^ Pa  is  the  mass  of  engulfed  air  at 

radius  R.  If  we  multiply  the  correction  term  by  the  cube  of  the 
radius  we  can  define  a new  radius  Z by 


Z = 


(26) 


where  M*  has  a definition  consistent  with  equation  (25) . We  define 
Z to  be  the  mass  corrected  radius. 

B.3  The  Waste  Heat  - Q 

Our  definition  for  the  energy  balance.  Equation  (18),  permits  us 
to  write  the  waste  heat  equation  as 


Q = Ae^  - /Pdv 


(27) 


since  the  total  energy  reduced  by  the  kinetic  energy  is  the  internal 
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energy.  We  can  use  Joule's  Law^^  to  write,  for  a perfect  gas, 


dE^  = C dT  . 
I V 


This  equation,  when  coupled  with  the  perfect  gas  law,  allows  us  to 
write 


Ae 


I 


AEi 

MASS 


P V 
a a 

Y-1 


(28) 


for  the  expansion  process. 


The  pres sure -volume  energy  is  found  using  GES: 


We  assume  that  our  expansion  process  is  adiabatic,  inviscid  and  non- 
conducting (i.e.,  isentropic) . This  is  permissible  since  we  are  not 
concerned  with  the  detailed  structure  of  the  shock  wave.  This  assump- 
tion sets  the  second  term  zero  since  the  assumption  results  in 
Pa  “ The  first  term  is 


P 


2 

a 

a a 


(i-y)y 


If  we  associated  the  final  state  with  Vg,  we  can  write  the  integral  as 
(dropping  the  subscript  I) 


P V 
a a 

1-Y 


Milne-Thomsorij  "Theoretical  Hydrodynamics,"  The  Mac  Mittan  Co., 
New  York,  1950 
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Again  making  use  of  the  isentropic  relation,  Equation  (23) , we  finally 
obtain  the  pressure-volume  energy 


W = /Pdv  = ^ 


(29) 


Combining  Equations  (27) -(29),  we  obtain  a dimensionless  waste  heat,  Q*: 


Q* 


-1 


(30) 


Lastly,  we  note  that  the  density  ratio,  P^/p  = D,  is  given  in 
terms  of  the  pressure  by  the  Hugoniot  relations  as 


Y+1 
P _ Y-1 

^a  p + 
r 


P +1 
r 


Y+1 

Y-1 


(31) 


We  are  interested  in  the  behavior  of  Q*  with  pressure.  For  small 
overpressures  Equation  (30)  can  be  expanded  by  using  the  binomial  theorem: 


(32) 


for  APj.  = (P  - Pa) Pa  Equation  (32)  states  that  in  the  acoustic 

wave  approximation  the  dimensionless  waste  heat  varies  as  (approximately) 
P^.  For  high  overpressures  (say  APj.  > 10)  Q*  is  found  to  more  approxi- 
mately follow  ® 


(22-L) (L-1) 

Q*  = 10  16  , L = Log^Q  (APp 


(33) 


and  gives  approximately  a linear  variation  between  Q*  and  AP^..  This 
is  illustrated  in  Figure  4. 
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Figure  4.  Variation  of  waste  heat  with  pressure 


B.4  The  Form  Factor 


We  introduce  one  more  concept  in  order  to  justify  the  QZQ  hypoth- 
esis: the  form  factor,  F.  To  quote  Porzel®:  "Probably  the  most 
important  single  experimental  fact  learned  about  explosions  in  the  past 
30  years  is  the  fact  that  they  scale,  and  over  enormous  ranges  of  yield, 
This  means  there  must  exist  a quantity  F which  is  not  unique  to  the  ex- 
plosion, requiring  a separate  calculation  for  each,  but  some  average 
energy,  a form  factor  common  to  all  explosions.  If  we  can  determine  it 
for  one,  we  can  determine  it  for  all  similar  explosions." 

In  general,  we  can  express  the  integrated  prompt  energy  as 


R 


Y = 4-n  f (W+K)r  dr  , 
0 

the  terms  of  which  can  be  arranged  as  follows: 


4 3 

Y = P 


r R 2 

I 3(W+K)  d(r/R)j 


or 


4 3 

Y = P 


f,  R ^ R ^ T 

^ / (W+K)r  dr^P  / r dr) 
o o 


(34) 


We  denote  the  bracketed  term  by  F,  the  form  factor.  It  is  a ratio  of 
the  average  energy  in  the  wave  to  the  peak  pressure.  The  integrated 
prompt  energy  can  thus  be  written 


Blast 

energy 


shock 

volume 


peak  pressure 
at  shock  wave 


F (35) 

average  energy  on 
interior  rel.  to 
peak  pressure. 


We  wish  to  determine  the  dependence  of  the  form  factor  on  the  radius, 
Since  Y ^ QR^  one  may  rewrite  Equation  (35)  as 


PF(P)  QR^/  Q 
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In  the  strong  shock  regime  P R"^  and  since  Q 'b  p r-3^  p is 
essentially  constant.  In  the  weak  shock  regime  the  spherical  wave  is 
very  nearly  planar  and  acoustic  theory  gives  P 'b  R“l.  But  Q 'b  p3 
in  the  weak  shock  regime  hence  F 'b  r-2. 

We  now  have  the  tools  necessary  to  justify  the  QZQ  hypothesis. 

B.5  The  QZQ  Hypothesis 


The  blast  energy  bal 
Y 

o 


ance  is 


R . 

47r  j (W+K)r  dr 
0 


+ 


R 

47r  j 
0 


Qr^dr. 


Initial  = prompt  energy  + waste  heat, 

yield 

Using  the  results  of  the  previous  section  Yq  can  be  written 


Y = 4^R^PF  + 4t:  / Qr^dr  . 
o 3 ^ 

o 

Differentiation  of  with  respect  to  R (noting  the  l/47r  d(Yo)/dr  = 0) 
gives 

QR^  + R^PF  + y R^  ^ (PF)  = 0 . 


Dividing  these  results  by  r2  P F we  have 


1 1 d 

3 RPF  dH 


(PF)  = 0 


or 


d(ln  PF) 
d(ln  R) 


In  the  strong  shock  regime  we  found  that  F was  essentially  constant 
and  Q 'b  p.  Therefore,  in  this  regime  Q/PF  'b  constant.  In  the  weak 
shock  regime,  Q 'b  p3^  p ^ r-1  and  with  F 'b  R“2  ratio  Q/PF  is  again 
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essentially  constant  (although  not  necessarily  the  same  constant 
as  that  in  the  strong  shock  regime) ; thus 


In  Q = In  constant  + In  PF. 


Denoting  3(1  + Q/PF)  by  q,  leads  to 

-q  = d[ln  Q ~ In  constant]/d  (In  R) 


or 


-q  = d(ln  Q)/d(ln  R)  . 

Integration  of  this  equation  produces  the  result  that 


QR^  = constant.  (^6) 

3 3 

Lastly,  we  note  that  R and  Z differ  by  an  additive  constant,  hence, 
their  derivatives  are  equal.  The  result  is  that  by  replacement  of 
derivatives,  integration  and  change  of  the  constant. 


constant. 


(37) 


We  require  estimates  for  the  constant  q in  the  strong  and  weak  shock 
regimes.  It  can  be  shown  that  the  above  equations  lead  to 


d(ln  Y)  _ d(ln  PF)  _ 
d(ln  R)  d(ln  R) 


Additionally,  Equation  (20b)  leads  to  the  suggestion  that 


d(ln  Y)  ^ 
d(ln  R) 


in  the  strong  and  weak  regimes.  In  the  weak  regime  all  of  the  prompt 
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energy  is  subject  to  waste;  however,  in  the  strong  regime  only  one 
half  the  prompt  energy  is  subject  to  waste.  Hence,  we  expect 


d(ln  Y) 
d(ln  R) 


r-1 


weak  regime, 
strong  regime. 


These  estimates  are  borne  out  by  detailed  calculations®  thus,  we  accept 
the  values  of  q as 


q = 


4.0  weak, 
3.5  strong. 
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c.  PEAK  OVERPRESSURE  CALCULATION  SCHEME  (POCS) 


The  necessary  technical  arguments  are  in  position  to  show  how 
to  make  a simple  scheme  for  calculation  of  spherical  blast  wave 
overpressure.  The  interior  analyses  give  rise  to  an  estimate  for 
the  initial  yield  and  an  initial  radius  for  the  charge  can  be  deter- 
mined by  1/6  of  the  propellant  mass  and  the  specific  gravity  of  the 
propellant.  All  that  remains  is  determination  of  the  constant  in 
the  QZQ  hypothesis. 

Since  QZ^  is  constant  for  all  Z we  shall  select  the  transition 
radius,  that  is,  the  radius  at  which  the  shock  wave  changes  from  strong 
to  weak,  as  the  point  at  which  the  constant  shall  be  evaluated.  Since 
the  initial  yield  and  radius  are  determined.  Equation  (20a)  can  be 
used  in  the  form 


Y 

o 


(39) 


where  Z^  is  the  transition  radius. 

The  constants  are  different  in  the  strong  and  weak  regimes  so  let 
us  write 


Strong: 

^1 

QZ  = A = Q^Z, 

Weak: 

QZ  ^ = B = Q^Z 

Equation  (39)  becomes  (using  dimensionless  energies  indicated  by  the 
asterisks) 


Y Z . 2-q  °° 

o ft  ^1  j f 

— A z dz  + 
i,  i 

Z Z^ 

o t 


4tt 


This  is  readily  integrated  to  give 


* 2-q 

B z dz  . 


3-q 


1 


(40) 
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since  q2  is  greater  than  3.  If  X is  defined  as  X = Z-^/Zq  then  we  can 
write  Equation  (40)  as 


* 


Y 

0 


3-q 

X 


1 


(41) 


This  equation,  once  is  specified,  can  be  solved  approximately  by  a 
variety  of  techniques.  (See,  for  example,  Schaumm's  "Numerical 

Analysis"^**)  , 


Equation  (30),  when  combined  with  Equation  (31),  gives  an  expres- 
sion for  Q*  in  terms  of  the  pressure  ratio,  P^,.  A unique  specification 
for  the  pressure  at  the  transition  point  is  not  possible  because  of  a 
variation  in  the  criteria  used  to  differentiate  strong  from  weak  shock 
waves.  For  example,  at  a pressure  ratio  of  3.8  x 10^  pascals  (3.8  atm.) 
the  sound  velocity  equals  the  material  velocity  behind  the  shock  wave 
and  this  is  a good  dividing  point  between  strong  and  weak  shocks.  Another 
reasonable  division  is  at  the  place  where  that  pressure  ratio  which 
separates  AP^  as  greater  or  less  than  Pa-  This  occurs  at  a pressure 
ratio  of  1 X 10^  pascal  (1  atm.).  Around  a pressure  ratio  of  2 x 10^ 
pascals  (2  atm.)  the  negative  phase  first  develops  thereby  preventing  any 
further  energy  from  propagating  from  the  interior  to  the  shock  wave. 

For  lack  of  any  definitive  criteria  for  transition  from  weak  to  strong 
shocks  we  select  a pressure  ratio  of  2 x 10^  pascals.  The  waste  heat 
at  the  transition  point  is  then  found  to  be 


Q^*  = .096. 


(42) 


In  conjunction  with  equation  (30),  (31),  and  (41),  the  pressure  ratio 
is  obtained  from 


.096(XZ  ) 
^ o 


^i 


(43) 


for  values  of  of  3.5  or  4 depending  upon  whether  Z is  greater  or  less 
than  Zt  (Z/Zo  ^ X). 


Soheid,  "Theory  and  Problems  of  Numeriaal  Analysis, " Sahawn ' s 
Outline  Series,  MoGraw-Eill  Book  Co.,  New  York,  (1968). 
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D.  ASPHERICAL  GEOMETRY 


The  previous  analyses  assume  spherical  symmetry  since  most  explo- 
sions can  be  idealized  as  emanating  from  a point  source.  The  gases 
exiting  from  a gun  would  probably  be  conical  or  cylindrical  in  shape 
owing  to  the  basic  shape  of  the  barrel  and  the  boundary  layer  buildup. 
We  wish  to  examine  the  effect  of  a cylindrical  charge  verses  a spheri- 
cal charge  upon  the  results.  We  do  this  by  replacing  the  idealized 
point  source  with  a line  source. 

The  basic  relation  for  prompt  energy.  Equation  (20a),  becomes,  for 
a line  source. 


Y = 2Trh  / Qzdz 

° z 

O 

where  quantities  are  assumed  constant  along  the  line  length,  h.  If  we 
assume,  for  the  moment,  that  the  exponents,  q,  are  the  same  in  both 
strong  and  weak  regimes  this  equation  integrates  to 


Y 

o 


Z-rrhKZ 
o 

q-2 


since,  as  before,  q is  greater  than  3.  Setting  K equal  to  the  values 
at  the  transition  radius  and  denoting  the  dimensionless  ratio 
(Z^/Z^)^  1 ^ ~ find  that,  for  equal  initial  yields 


or 


ZTThQ^Zo^ 


4TrQ^Zc 


q-2 


q-3 


2 Zq  ((i-2) 


Since  the  initial  volumes  must  be  equal,  i.e.. 


one  has 


2 


sVq-Z) 

4;(q-3)  ’ 
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or 


3 £-_2 
2 q-3 


1/q 


We  are  left  with  estimating  the  initial  radius  of  the  clyinder 
with  respect  to  that  of  the  sphere.  In  the  next  chapter  we  will  show 
that  the  ordered  motion  is  approximately  equal  to  the  random  motion, 
hence,  energy  transport  along  the  axis  of  the  barrel  would  be  greater 
than  the  transverse  flux  by  a factor  of  two.  This  suggests  an  aspect 
ratio,  h/Z,  of  2 and,  since  the  volumes  must  be  equal,  we  obtain 


I 2, 


or 


Zo./Zos 


.8735. 


Consequently, 


.8735 


This  function  is  plotted  in  Figure  5 and  shows  that  the  transition 
radius  is  from  1.5  to  2.5  times  farther  away  when  the  initial  energy 
release  is  cylindrical  rather  than  spherical  in  shape. 


1 £li- 

2 q-3 


i/q 


(44) 
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q 


Figure  5.  Cylindrical  fcharge  contribution  to  transition  radius  scaling 
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li. 


i’ULSb  LliNGTll  DHTLiRM  I NATION 


The  traditional  analysis  to  determine  the  blast  pulse  length  uses 
an  assumed  pressure  pulse  time  history  and  empirical  blast  data.  Pulse 
shapes  ranging  from  polynomials  (including  instant-rise/linear-decay 
triangles)  to  the  product  of  polynomials  and  exponentials  are  commonly 
used.  A classical  example  is  the  Friedlander  equation  as  found  in 
Baker® : 


P(t) 


P.  + 
a 


PJ  1 


exp(-bt/t_^). 


In  the  usual  analysis  the  four  constants  Pa,  Pp,  t+  and  b are  determined 
(as  a function  of  charge  weight  and  distance  from  the  charge)  empiri- 
cally. 

Another  approach  is  that  of  Theilheimer® . He  uses  a pulse  shape 
of  the  form 


P(t)  = P + P„  exp(-t/e  ) 
a p 


where  the  time  constant  is  defined  by 


0 


dP/dtJ 


evaluated  at  t = 0+.  Using  this  pulse  shape  and  the  partial  differen- 
tial equations  that;  govern  the  conservation  of  mass  and  momentum; 
define  the  speed  of  sound;  and  describe  the  change  of  (energy)  state 
in  an  adiabatic  expansion,  Theilheimer  finds  the  partial  derivative  of 
pressure  with  respect  to  time  as 


IP 

3t 


2pua 


(U-u)  . § 


a + u(U-u) 


du  2 

* dis  “ 


(U-u)' 


The  UTE  could  be  used  to  determine  gradients  of  u and  P with  respect 
to  R and  then  to  3P/3t  and  from  it,  the  time  constant  for  decay  of  pres- 
sure at  fixed  distances  from  the  blast  center.  However,  the  data  cur- 
rently available  on  durations  does  not  justify  this  close  examination. 

In  order  to  apply  the  UTE  to  pulse  length  calculation  without  resort 
to  a complex  computation  (which  is  not  justified  by  the  quality  of  the 
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data)  we  note  that  the  theory  as  presented  is  not  time  dependent; 
however,  it  can  be  extended  for  estimating  the  pulse  length.  We  begin 
by  observing  that  a time  can  be  determined  using  a characteristic 
length  associated  with  the  energy  volume  and  the  material  velocity. 
Additionally,  one  can  relate  the  pressure  volume  energy  to  the  prompt 
energy. 

Let  us  define  the  pulse  length  by 


At  = 5/u 


where  ^ is  representative  of  the  volume  occupied  by  the  pressure  energy. 
We  define  5 by 


4 Av  ^3 

T TT  R . 

3 V 

o 


In  general  then  (from  Equation  (35)) 


Yh(e) 


4 0 ,3_ 

— TT  — 5 P 
3 Av 


where  h(e)  modifies  the  yield  to  include  only  that  portion  which  appears 
as  pressure  energy. 

The  volume  ratio  (because  of  the  per  unit  mass  definition)  becomes 


Av  ^ V D-1 
V 

o 

where  D = p/po-  The  form  adopted  for  h(e)  is  arrived  at  as  follows: 

The  prompt  energy  is  defined  as  the  static  overpressure  energy 
(internal  energy,  ej)  plus  the  kinetic  energy.  Since  the  kinetic 
energy  is  not  detected  by  a static  pressure  probe  we  must  reduce  the 
prompt  energy  by  the  ratio: 


Ae, 


AOj  + A(u  /2) 
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This  can  be  written  as 


1 

^ ^ A(uV2) 

Aei 

and,  since  we  have  assumed  a perfect  gas,  constant  total  enthalpy 
allows  us  to  write  the  ratio  as 


1 


The  ratio  of  the  enthalpy  to  internal  energy  is  recognized  as  the  ratio 
of  the  specific  heats,  yj  hence,  the  correction  factor  for  the  prompt 
energy  becomes 


h(e) 


1 

1+y  * 


The  pulse  length  becomes 


= i CD-1)Y 

^ j iT(l+y)PD 


(46) 
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Ill  - COMPARISON  OF  THEORY  AND  EXPERIMENT 


The  theory  which  we  have  presented  raises  many  questions  concerning 
the  assumptions  and  the  ultimate  validity  of  such  an  approach.  In  this 
section  we  shall  examine  data  for  three  guns  in  the  Army  inventory  and 
compare  these  data  to  our  theory.  The  guns  considered  and  the  relevant 
characteristics  are  given  in  Table  I.  (We  apologize  for  the  inconsis- 
tent units;  hov/ever,  this  is  the  way  data  is  usually  presented^^j . 


Table  I 


(a)  Gun 


Designation 

Type 

Barrel 
Length (m) 

Dia . (mm) 

Groove 
Height (m) 

Propellant 

XM204 

towed 

howitzer 

3.55 

105 

3.5x10-2 

M30AI 

MII0E2 

self-propelled 

howitzer 

6.93 

203 

3.66x10'^ 

MI88E2 

MI  07 

self-propelled 

8.95 

175 

3.66E-3 

M6 

artillery  gun 
(b)  Propellant 


Specific 

Gravity(g/cc) 

Specific 

Energy(cal/g) 

Chamber 

Temp(K) 

Chamber 
Press . (psi) 

Weight/ zone 
(lbs) / (non-dim) 

M30AI 

1.66 

975 

3040/2450 

54000 

4.42/8 

MI88E2 

1.66 

975 

3040/2450 

31000 

38/8 

M6 

1.58 

758 

3040/2450 

46000 

57/3 

Reiahardj  A.R.  Downs,  "A  Compendium  of  Field  Avtillevy  Facts  - 
Organization,  Tactics,  Operations,  Weapon  Systems  and  Terminology , " 
BRL  Report  Do.  1759,  Feb  1975.  (AD  #8002431) 

"Interior  Ballistics  of  Guns, " Engineering  Design  Handbook,  AMO 
Pamphlet  Ho.  AMCP  706-150,  Feb.,  1965. 
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Table  1 (Cont) 


(c)  Projectile 

Gun 

Weight (lbs) 

Muzzle  Vel (ft/sec) 

XM204 

33 

2133 

M110E2 

200 

2330 

Ml  07 

147.5 

3000 
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A.  APPLICATION  OF  SPHERICAL  EXPLOSION  THEORY  TO  MUZZLE  BLAST 
OVERPRESSURE  PREDICTION 


The  overpressure  theory  of  the  previous  chapter  is  not  directly  ap- 
plicable to  the  muzzle  blast  problem.  In  this  section  we  shall  develop 
those  adjustments,  based  upon  preliminary  examination  of  the  data,  to 
make  the  principal  result,  the  QZQ  hypothesis,  applicable.  In  this 
section  we  shall  address  three  areas:  (1)  the  effect  of  the  ordered 
motion  of  the  propellant  gases,  (2)  revised  estimates  for  the  expo- 
nents in  the  QZQ  hypothesis  and  (3)  the  effect  of  the  cylindrical 
vs.  spherical  initial  shape.  Finally,  we  shall  illustrate  our  findings 
by  using  the  175mm  Ml 07  artillery  gun. 

A.l  The  Moving  Charge  Effect 

Energy  transport  in  any  direction  is  determined  by  the  sum  of 
ordered  and  random  motion,  each  contributing  according  to  their  respec- 
tive velocities.  The  gases  issuing  from  a gun  have  velocities  given  by 
the  projectile  muzzle  velocity,  which  is  on  the  order  of  3 times  the 
speed  of  sound;  hence,  it  is  instructive  to  examine  the  relative  signi- 
ficance of  these  motions.  Experimental  investigations  into  the  effect 
of  motion  on  blast  overpressures  were  performed  by  Patterson  and  Wenig^^, 
Armendt^®  and  Armendt  and  Sperrazza^^,  Their  results  showed  that  blast 
overpressures  were  measureably  greater  at  a given  distance  from  detona- 
tion in  the  direction  of  motion  as  compared  to  the  transverse  direction. 
We  expect  a similar  effect  in  muzzle  blast  and  in  the  following  we  give 
a simple  procedure  for  including  this  effect. 

We  begin  with  a comparison  between  the  energies  of  random  and 
ordered  motion.  The  kinetic  energy  of  ordered  motion  is  simply 


1 2 

2 muzzle 


Patterson^  J.Wenig^  Blast  Measurements  Around  Moving 

Explosive  Charges^^^  BRL  Memorandum  Report  Eo.  76? ^ Mar.^  1954,  (AD  #33173) 

18 

B.F.  Armendt 3 "Air  Blast  Measurements  Around  Moving  Explosive 
Charges,  Part  II,"  BEL  Memorandum  Report  No.  900,  May,  1955.  (AD  #71277) 

19 

B.F.  Armendt,  J.  Sperrazza,  "Air  Blast  Measurements  Around  Moving 
Explosive  Cha.rges,  Part  III,"  BRL  Memorandum  Report  No.  1019, 

July,  1956  CAD  #114950) 
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which,  for  the  Ml 07,  is 


j/kg 


or 


ordered 


915  m/sec  . 


The  energy  in  random  motion  is  given  by 


T 


per  degree  of  freedom  where  k is  Boltzmann *s  constant.  Table  1 gives 
the  temperatures  as:  isochoric  - 3040  ®K,  isobaric  - 2450  °K.  We 
expect  the  actual  temperature  to  be  somewhere  between  these,  hence, 
for  convenience,  let  us  assume  the  temperature  to  be  2750^K  and  an 
average  molecular  weight  of  28.  The  energy  in  random  motion  becomes 


i k T = 4.1  X 10^  j/kg 


which  gives  a random  velocity  of  approximately  904  m/sec. 

These  results  show  that  energy  transport  along  the  axis  of  the 
barrel  is  approximately  twice  that  in  the  transverse  direction.  We 
interpret  this  as  meaning  that  pressure  probes,  one  at  a distance  L 
perpendicular  to  the  barrel  axis  will  sense  the  same  pressure  as  one 
aligned  along  the  axis  at  a distance  L/2.  For  simplicity,  we  adjust 
all  lengths  by  the  factor 


1 

1 + cos  0 


where  6 is  the  angle  between  a unit  vector  along  the  barrel  axis  (at 
the  muzzle)  and  a position  vector  from  the  muzzle  to  the  pressure 
probe.  This  function  has  the  characteristic  of  varing  smoothly  from 
1/2  along  the  barrel  axis  to  unity  normal  to  the  axis. 
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Figures  6a  and  6b  show  the  effect  of  this  modification.  F’igure 
6a  shows  the  data  as  measured.  The  symbols  denote  gun  elevations 
of  : + 100  mils,  □ 620  mils  and  X2lll25  mils  (military  mils  where  6400 
mils  equal  360^) . The  variation  in  each  symbol  is  a result  of  differ- 
ent azimuth  angles  (0®,  30®,  60®,  90®).  Application  of  relation  (47) 
clearly  results  in  a more  reasonable  correlation  of  overpressure  with 
distance  from  the  muzzle.  For  this  reason  all  data  is  adjusted  accord- 
ing to  this  prescription. 

A. 2 Revised  Estimates  for  q 

In  Section  II. B. 5 values  of  3.5  and  4 were  obtained  for  the  ex- 
ponents, q,  in  the  strong  and  weak  shock  regimes  respectively.  These 
values  were  predicted  upon  the  scaling  laws  which  state  that  P ^ R"1 
in  the  weak  regime. 

The  previous  section,  III.A.l,  shows  that  there  is  an  effect 
due  to  the  forward  velocity  of  the  gases.  Since  the  q’s  relate  the 
conversion  of  prompt  energy  to  waste  heat,  and  since  the  ordered  motion 
is  approximately  equal  to  the  random  motion,  we  propose  that  only  half 
of  the  pressure  volume  energy  is  available  for  conversion  to  waste  heat. 
The  consequence  is  that  in  the  strong  regime  q 3.25.  Additionally, 
since  the  muzzle  overpressures  are  relatively  low  (compared  to  ii.E. 
explosives)  the  transition  radius  is  small  compared  to  the  distances 
of  interest.  Hence,  the  motion  effect  has  not  disappeared  and  we 
plausibly  expect  the  weak  regime  q to  be  approximately  3.25  also. 

This  conclusion  is  demonstrated  in  Figure  7 where  we  have  plotted 
-In  Q vs.  In  R (or  In  Z since  they  are  approximately  equal  far  from  the 
muzzle)  for  the  adjusted  data  of  Figure  6b.  The  straight  line  is  a 
least  squares  regression  fit  of  the  form  ax^.  The  fit  is  to  the 
increasing  data;  the  points  at  the  greatest  distances  were  not  included 
since  they  are  at  essentially  constant  pressure  and  do  not  reflect  the 
dynamics  of  the  processes.  For  these  data  a slope  of  approximately  3.2 
was  found.  This  is  consistent  with  our  assertion  and  demonstrates  that 
the  division  of  energy  in  crossing  the  shock  wave  is  significantly 
affected  by  the  existence  of  the  ordered  motion. 

In  our  analyses  we  shall  use  equal  q's  of  value  3.25. 

A. 3 Initial  Cylindrical  Shape 

In  Section  D of  the  previous  chapter  we  developed 


47 


OVERPRESSURE  RATIO  - AP.  OVERPRESSURE  RATIO 


DATA  ADJUSTMENT  FOR  MUZZLE  VELOCITY  EFFECT 
M107  ARTILLERY  GUN  (175MM) 


DISTANCE  FROM  MUZZLE  - METRES 
FIGURE  6 
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VARIATION  OF  WASTE  HEAT  WITH  DISTANCE 
M107  ARTILLERY  GUN  (175MM) 
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a modification  to  the  theory  because  of  a proposed  initial  cylindrical 
shape  to  the  charge.  For  a q of  3.25  we  find  that 


cyl 


1.624  (ZJ 

sph 


When  the  transition  radius  is  near  the  muzzle  we  still  expect  the 
cylindrical  shape  to  be  present  since  sufficient  time  has  not  elapsed 
for  these  asymmetrical  effects  to  be  washed  out.*  That  the  cylindrical 
transition  radius  is  more  appropriate  is  shown  in  Figure  8. 


\ 


*For  an  exit  Mach  number  of  2.3  and  transition  radius  of  1 metre  the 
elapsed  time  is  of  the  order  of  1 msec. 
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DISTANCE  FROM  MUZZLE  - METRES 


B.  COMPARISON  OF  THEORY  AND  EXPERIMENT 


In  this  section  we  will  apply  our  theory  to  the  guns  described 
at  the  beginning  of  this  section.  In  each  case  we  have  also  pro- 
vided a least  squares  regression  fit  of  the  form  y = a to  the 

adjusted  data.  The  data  for  the  M110E2  and  the  Ml 07  were  obtained 
from  the  Material  Test  Directorate  at  the  Aberdeen  Proving  Ground, 

Md.^^  The  XM204  data  were  obtained  from  Westline *s  report^.  We 
note  that  the  APG  data  error  is  claimed  to  be  a maximum  of  5%.  We 
also  note  that  for  clarity  in  certain  instances,  data  which  did  not 
add  further  information  were  not  plotted. 

B.l  Muzzle  Blast  Overpressure 

Figure  8 in  the  previous  section  shows  our  prediction  for  the 
175mm  M107  artillery  gun.  For  the  adjusted  data  the  regression  fit 
yielded  an  exponent,  b,  of  -1.29.  The  measured  data  are:  + - 100  mils 
quadrant  elevation  (Q.E.),n-  620  mils  Q.E.,  and  - 1125  mils  Q.E. 
Additionally,  each  symbol  includes  variation  in  the  azimuth  at  values 
of  0®,  30®,  60®,  and  90®.  Two  probes  were  located  at  each  of  the 
angles  on  circular  arcs  of  radii  9.15  and  12.2  metres  from  the  muzzle 
when  at  a Q.E.  of  100  mils. 

Within  the  data  field  a maximum  disagreement  (in  overpressure  ratio) 
between  the  theory  and  data  is  found  to  be  approximately  +15%.  This 
corresponds  to  a 3.9%  error  in  the  absolute  pressure. 

Figures  9a  and  b show  the  data  and  theoretical  prediction  for  the 
M110E2  self-propelled  howitzer.  Figure  9a  shows  the  measured  data  with 
the  symbols  denoting  the  same  conditions  as  with  the  Ml 07. 

Figure  9b  shows  the  adjusted  data,  the  theoretical  prediction  and 
the  least  squares  regression  fit.  The  regression  fit  gave  an  exponent, 
b,  of  approximately  -1.26.  In  this  fit  the  two  flagged  points  were  not 
included.  These  points  correspond  to  a Q.E.  of  100  mils  with  the  probes 
located  directly  in  front  of  the  muzzle  (0®  azimuth).  These  abnormally 
low  overpressures  may  be  attributed  to  the  losses  associated  with  the 
flow  patterns,  i.e.,  the  normal  shocks  associated  with  the  Mach  disks, 
and  hence,  are  not  characteristic  of  the  overall  decrease  of  pressure 
with  distance.  This  assertion  cannot  be  validated  since  similar  data 
were  not  observed  with  the  other  guns  examined.  We  can  state  that 
since  the  basic  reason  for  this  paper  is  predication  of  overpressures 
for  safety  reasons,  these  low  values  are  of  no  consequence  for  the  work 
at  hand. 


Laoijy  Private  Corrmtni cation. 
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OVERPRESSURE  RATIO  - AP.  OVERPRESSURE  RATIO 


MUZZLE  BLAST  OVERPRESSURE 
Ml  10E2  SELF-PROPELLED  HOWITZER  (203MM) 


FIGURE  9 
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For  these  M110E2  data  the  theoretical  prediction  is  upwards  of 
+25%  in  error  for  the  overpressure  ratio  (compared  to  the  regression 
fit).  This  corresponds  to  approximately  5.2%  in  absolute  pressure. 

Figures  10  a and  b show  the  data  for  the  XM204  towed  howitzer. 

Figure  10a  shows  the  measured  data.  These  data  were  all  obtained  at 
a quadrant  elevation  of  26.7  mils  and  the  data  represent  variation  in 
the  azimuth. 

Figure  10b  shows  the  adjusted  data.  For  these  data  the  theore- 
tical prediction  and  the  regression  fit  are  basically  indistinguish- 
able. This  figure  also  nicely  illustrates  how  the  data  collapse 
about  the  90®  azimuth  value. 

B.2  Muzzle  Blast  Pulse  Length 

The  data  in  this  section  represent  the  positive  phase  duration. 

That  is,  the  measured  data  represent  the  time  duration  between  the 
start  of  the  pulse  and  the  point  where  the  initial  pulse  returns  to 
ambient  pressure.  The  range  component  of  the  duration  data  has  been 
adjusted  in  the  same  way. (and  for  the  same  reason)  as  the  range  com- 
ponent of  the  pressure  data.  But,  unlike  the  pressure  data  the  pulse 
length  is  adjusted  in  addition  to  the  range.  We  argue  that  an  observer 
standing  in  front  of  the  muzzle  will  see,  in  time,  a shorter  pulse 
because  of  its  velocity  than  an  observer  standing  at  90®  to  the  muzzle 
unit  normal.  We  have  provided  the  regression  fit  for  these  data; 
however,  there  is  sufficient  scatter  in  the  data  that  both  authors 
and  readers  are  wise  to  refrain  from  drawing  any  comparisons  between 
the  theory  and  the  fit. 

Figures  11  a and  b show  the  pulse  length  for  the  M110E2  and 
Figures  12  a and  b show  the  data  for  the  Ml 07.  The  symbols  correspond 
to  the  same  conditions  as  noted  for  the  overpressure  data.  Our  theo- 
retical predictions  generally  fall  along  the  top  of  the  measured  data. 

We  have  not  examined  the  XM204  pulse  length  data  because  of 
difficulties  in  interpreting  the  exact  nature  of  the  data. 

Lastly,  Figure  13  shows  the  M107  overpressure  data  again,  this  time 
including  bars  to  denote  the  stated  error  of  +5%.  We  show  this  to 
reaffirm  the  discussion  which  compared  our  theory  to  the  regression  fit. 
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OVERPRESSURE  RATIO  - AP,  OVERPRESSURE  RATIO 


MUZZLE  BLAST  OVERPRESSURE 
XM204  TOWED  HOWITZER  (105MM) 


DISTANCE  FROM  MUZZLE  - METRES 
FIGURE  10 
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MUZZLE  BLAST  PULSE  LENGTH 
M110E2  SELF-PROPELLED  HOWITZER  (203MM) 


b.  ADJUSTED  DATA 
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DISTANCE  FROM  MUZZLE  - METRES 
FIGURE  11 
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PULSE  LENGTH  - MSEC  PULSE  LENGTH  - MSEC 


MUZZLE  BLAST  PULSE  LENGTH 
M107  ARTILLERY  GUN  (175MM) 
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b.  ADJUSTED  DATA 


DISTANCE  FROM  MUZZLE  - METRES 
FIGURE  12 
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C.  DISCUSSION 


In  view  of  our  results,  specifically  with  regard  to  the  modification 
for  cylindrical  geometry,  we  can  form  a tentative  picture  of  the  gross 
features  of  the  muzzle  blast  flow  field. 

Perhaps  the  most  unrealistic  characteristic  of  the  assumed  cylin^ 
drical  shape  is  the  requirement  of  constant  pressure  along  the  length 
of  the  cylinder.  An  examination  of  the  expected  processes  in  the  flow 
field  show  that  this  may  well  be  true.  Previously  we  spoke  of  the 
inter-relationship  between  the  static  and  dynamic  pressures  dla 
Bemoull'i.  In  those  discussions  we  used  the  Bernoulli  principle  to 
arrive  at  the  projectile  base  static  pressure.  Once  the  gases  have 
left  the  muzzle  the  reverse  process  must  occur.  That  is,  the  high 
dynamic  pressure  of  ordered  motion  must  revert  to  static  pressure  as 
the  gases  decelerate  to  zero  velocity.  This  velocity  decrease  will 
occur  along  the  axis  of  the  gun  hence  we  should  expect  a production  of 
pressure  energy  along  this  line- -forming  a more  cylindrical,  rather  than 
spherical,  shape. 

This  energy  conversion  produces  a second,  equally  important  effect. 
Classical  explosion  theory  assumes  an  instantaneous  point  source  explo- 
sion. Porzel's  theory  admits  of  afterburning*  by  permitting  q's  below 
his  suggested  values  but  we  found  that  even  reducing  the  exponent’s 
value  to  the  minimum  (lower  bound  of  3)  did  not  produce  the  necessary 
change  in  the  slope  of  the  pressure  curve.  This  production  of  pressure 
energy  persists  for  such  sufficiently  long  times  compared  to  ordinary  ex- 
plosions, that  this  pressure  energy  production  continues  to  feed  the 
shock  wave  long  after  the  weak  shock  regime  has  been  entered,  as  is 
demonstrated  in  Figure  14. 

Figure  14  shows  the  absolute  peak  pressure  plotted  against  the 
range.  We  applied  a least  squares  linear  regression  fit  to  the  de- 
creasing pressure  data  and  found  a slope  of  approximately  -.5.  This 
is  half  that  determined  for  standard  spherical  explosions  and  is 
indicative  of  some  more  or  less  continuous  source  of  pressure  energy 
addition  to  the  shock  front. 

One  might  alternatively  view  this  in  light  of  the  concepts  of  prompt 
energy  and  waste  heat.  We  have  noted  that  the  kinetic  energy  is  not 
subject  to  waste.  In  the  muzzle  blast  problem  we  have  a large  reser- 
voir of  kinetic  energy  due  to  the  initial  high  velocity  of  the  gases; 
because  of  this  resevoir  the  energy  division  suggested  by  Rankine- 
Hugoniot  relations  is  not  valid. 


*A  process  of  continued  energy  addition  to  the  shock  front. 


LINEAR  REGRESSION  FIT  TO  MUZZLE  BLAST  PEAK  PRESSURE 
M107  ARTILLERY  GUN  (175MM) 
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DISTANCE  FROM  MUZZLE 


D.  SUMMARY  CONCLUSIONS 


In  general  the  a priori  muzzle  blast  theory  is  viable.  Overpressure 
prediction  agrees  with  the  adjusted  data  to  within  the  stated  error  in 
data.  The  pulse  length  prediction  generally  follows  the  upper  (longer) 
pulse  length  data  which,  for  safety  requirements,  is  a most  desirable 
attribute.  We  note  that  there  is  sufficient  scatter  in  the  pulse  length 
data  that  further,  more  refined  experimentation  is  warranted  before  com- 
plete confidence  is  obtained.  In  this  regard,  more  detailed  study  of 
Porzel*s  theory  in  deriving  the  constant  may  be  justified,  as  suggested 
in  Section  II. E. 

At  the  outset,  we  suggested  that  a detailed  knowledge  of  the  flow 
field  is  not  necessary.  This  has  been  demonstrated  for  the  data 
examined.  We  recognize  that  all  the  data  occurred  in  the  weak  shock 
regime;  had  strong  shock  regime  data  been  available,  we  would  have 
expected  a departure  of  this  theory  from  experiment  - a consequence 
of  the  detailed  flow  field.  We  would  expect  the  theory  to  overpredict 
since  there  are  losses  in  the  muzzle  flow  field  not  accounted  for  in 
the  UTE.  We  do  note  that  the  transition  distance  separating  the 
strong  from  the  weak  regime  was  consistently  less  than  one  metre, 
a data  source  for  R < Im  would  not  likely  be  found. 

We  have  examined  three  guns  of  sufficient  variation  in  length, 
loading  etc. , to  span  the  data  field  for  large  guns  and  therefore, 
the  theory  can  be  reasonably  applied  to  any  large  gun  in  the  Army 
inventory,  but  only  within  the  constraints  of  our  analysis. 

We  have  not  examined  the  small  gun  effects  - an  area  we  leave  for 
further  application  of  the  theory.* 

The  most  striking  result  of  this  investigation  is  the  difference 
in  the  behavior  of  the  shock  overpressure  with  distance  between  the 
guns  examined  and  that  which  results  from  conventional  point  source 
explosions.  The  energy  partition  concepts  suggest  that  a significant 
reduction  in  overpressure  can  be  achieved  by  introducing  some  effective 
means  for  converting  this  high  dynamic  pressure  to  static,  or  random 
motion,  early  in  the  propagation  phase.**  The  implementation  of  such 
a device  would  probably  increase  the  overpressure  near  the  muzzle  but, 
since  higher  overpressures  mean  higher  losses,  less  energy  would  be 
available  in  the  shock  as  it  leaves  the  immediate  vicinity  of  the 
muzzle.  Considerable  reduction  in  the  overpressure  can  be  achieved  if 
the  range  dependence  can  be  made  to  more  closely  approximate  the  point 
source  dependence,  i.e.,  R"l. 


*See  next  section  for  recommended  future  studies. 

**In  small  arms,  noise  suppressor  (silencers)  do  exactly  this  through 
a throttling  process. 
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This  report  presents  a theory  which  is  applicable  to  the  muzzle 
blast  problem  for  the  situations  considered.  Several  variations  were 
not  examined.  The  next  section  suggests  future  studies  with,  in  some 
instances,  possible  techniques  for  including  the  additional  effects. 


/ 
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E.  SUGGESTED  FUTURE  STUDIES  - AREAS  NOT  CONSIDERED  IN  THIS  STUDY 


1.  Small  Arms  - We  know  of  no  a priori  reason  why  the  theory 
would  not  be  directly  applicable.  One  possible  exception  would  occur 
if  the  gun  L/D  becomes  less  than  the  L/D  for  the  boundary  layer  choke. 
Should  this  occur,  we  believe  the  muzzle  velocity  of  the  bullet,  and 
hence  the  exhaust  gases,  will  be  indicative  of  the  energy  reduction 
from  the  chamber  pressure. 

2.  Muzzle  Brakes  - We  view  this  as  strictly  an  energy,  mass,  and 
flow  symmetry  diverter.  It  is  probable  that  most  of  this  effect  can  be 
included  in  MEZ  and  the  initial  yield  estimate.  However,  an  appropriate 
trigonometric  function  describing  the  additional  variations  in  angle 
about  the  barrel  axis  will  be  required, 

3.  Reflected  Shocks  - (a)  from  the  ground.  One  possibility  is  to 
double  the  initial  yield  input  to  the  calculation.  Another  is  to 
provide  an  accelerated  rate  of  conversion  of  kinetic  energy  of  ordered 
motion  to  random  motion  at  the  various  ranges,  (b)  from  structures 
such  as  gun  mounts.  Detailed  analysis  near  the  structure  would  be 
very  complicated  because  of  the  geometry  of  the  structure  itself.  Away 
from  the  structure,  however,  local  spatial  variations  should  smooth 
out.  Perhaps  introducing  a cylindrical  explosion  of  shape  equivalent 
to  the  structure  would  provide  a start.  As  with  the  current  problem, 
an  estimate  for  the  equivalent  yield  would  probably  be  the  more  diffi- 
cult task. 

4.  Lastly,  we  would  like  to  see  more  refined  tests  so  a better 
evaluation  of  this  theory  can  be  made. 
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Appendix  - Computational  tools 

At  the  outset  we  sought  a theory  which  was  simple  to  use.  In 
this  appendix  we  supply  two  different  schemes  for  calculating  the 
muzzle  blast  overpressure  and  pulse  length.  We  are  presenting  these 
different  schemes  because,  depending  upon  the  user  and  the  situation, 
diff®i'6nt  levels  of  versatility  may  be  required.  These  two  schemes 
and  their  scope  are: 

(1)  A code  for  the  Texas  Instruments  SR-52  programmable  pocket 
calculator.  In  this  version  all  assumptions,  i.e.,  equal  q,  cylindri- 
cal symmetry  etc.  are  incorporated  and  combined  into  "lumped" 
constants.  Changing  any  of  the  values  requires  considerable  reprogram- 
ming.  To  conserve  space  in  the  machine  we  have  applied  "fits"  to  the 
transcendental  equations  and,  additionally,  the  pressure  ratio  is 
input  and  the  radius  determined  rather  than  vice  versa. 

(2)  A code  written  in  BASIC  suitable  for  use  on  a minicomputer. 
Complete  flexibility  is  permitted  in  this  code.  This  program  was 
developed  on  a NOVA  minicomputer  with  8k  memory  (of  which  the  BASIC 
translator  occupies  approximately  4k) . 

Also,  a series  of  nomographs  will  be  provided  under  separate  cover 
where  quantities  are  considered  together  as  combined  variables.  Their 
numerical  values  are  left  to  the  user.  In  the  nomographs  we  have  assumed 
the  exponents,  q,  to  be  equal  and  also  that  the  user  has  available 
some  sort  of  calculating  machine  whether  it  be  a pocket  calculator  or 
a slide  rule. 

These  computational  schemes  are  designed  to  be  used  with  the  metric 
system  of  units,  specifically  the  mks  system.  To  facilitate  use  of 
these  codes  by  those  not  accustomed  to  working  in  these  units  we  pro- 
vide the  following  table: 


Table  App-1 


Units  Given 

Times 

Metric  Units 

Length  ft. 

3.048x10-1 

metres  - m 

Wgt/mass  lbs. 

4.536x10-1 

Kilograms  - kg 

Pressure  psi 

6. 894x1 0^ 

Pascals  (Pa)  - 
J/m^ 

Energy  ft. -lb. 

1.356 

Joules  - J 

Velocity  ft. /sec 

3.048x10-1 

metres/sec  - 
m/sec 
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App-1  Code  for  the  Texas  Instruments  SR-52  Pocket  Calculator. 

This  code  is  contained  on  two  cards.  The  first  program  accepts 
the  input  data  and  calculates  the  yield  and  radius  of  the  charge.  The 
second  program  accepts  the  pressure  ratio  and  determines  the  distance 
from  the  muzzle  at  which  the  pressure  ratio  occurs  and  the  pulse  length 
at  that  point.  The  code  makes  use  of  fits  to  the  complicated  functions 
to  permit  an  estimate  of  the  correct  answer. 

The  analysis  is  as  follows:  (Quantities  followed  by  (1)  are  input) 

(1)  The  kinetic  energy  of  the  projectile  is  obtained  from  the  mass  (1) 
and  velocity  (1)  of  the  projectile. 

(2)  The  total  energy  is  obtained  from  the  mass  (1)  and  specific 
energy  (1)  of  the  propellant. 

(3)  . The  difference  between  (2)  and  (1)  becomes  the  maximum  available 
energy. 

(4)  The  chamber  overpressure  ratio  (I)  is  used  to  determine  the  choke 
overpressure  ratio  from 


(AP  ) 

Choke 


exp 


.92  ln(AP  ) -2.14 

Chamber 


(5)  The  barrel  length  (I),  barrel  diameter  (I)  and  the  groove  height 
(1)  are  used  to  calculate  the  rough  wall  tube  loss  from 


If  (L/D)(.j^Q]^g  is  greater  than  (L/D)gun  the  program  halts  with  a blink- 
ing display  - this  analysis  does  not  include  this  situation. 

(6)  The  program  then  determines  the  yield  via 

Yo  = fCAP^)  / 

/ Chamber 


6 X (AP  ) 
r 


Muzzle 


"Avail 


-1 
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(7)  The  mass  correction  for  MEZ  (Section  II. B. 2)  is  determined  by 


M'  = 8.2  X 10" V -m^ 

Prop 


(8)  The  specific  gravity  of  the  propellant  (I)  is  used  to  calculate 
the  inital  radius  from 


R ^ 
o 


25p  Prop 
P 


1/3 


-m  . 


(9)  The  mass  corrected  initial  radius  is  obtained  from 


'o  = [“o'  - «’] 


1/3 


-m  . 


(10)  The  transition  radius,  Z^,  is  found  from 

7 Y 1/3-25 

Z = Z Z f _m 

t o[Z^J  oL 


O'  '16^Qt^o‘ 


-2 


where  the  value  of  Q^,  9.55x10  , is  found  assuming  an  overpressure 

ratio  of  2 at  the  transition  point. 

This  concludes  the  calculations  on  card  1.  Card  2 calculates 
the  following; 

(11)  The  uncorrected  transition  radius,  R^,  is  obtained  from 

r 3 “I  1/3 

\ ■ “'J 

(12)  The  constant  in  the  QZQ  hypothesis  is  obtained  from 


- Qt\ 


3.25 


for  a of 


9.55x10"^. 
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(13)  The  program  now  accepts  a pressure  ratio  (I)  and  calculates 
the  distance  from  the  muzzle  at  which  this  pressure  ratio  occurs  via 


and 


6P  +1 


R = 


Z^-  M' 


1/3 


X 1.624 


(14)  The  program  then  calculates  the  pulse  length  from 


u = 280 


(P^-l) 


/D-1 

I D 


1 1/2 


At  = - 
u 


u 


m/sec 

1/3 


|lT(P)(l+Y)  (^)J 


5A(D-1)Z 


-.25 


D(P) 


1/3 
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The  user 

instructions,  required  inputs  and  displayed 

quantities  follow; 

Step 

Procedure 

Enter  Press 

Display 

Card  1 

1 

Load  side  A 

Clr  2nd 

Read 

2 

Load  side  B 

2nd 

Read 

3 

Initialize 

Clr  A 

0. 

4 

Kinetic  En.  of 

Proj . mass  - kg 

RUN 

M/2 

Pro  j . 

5 

Proj . Vel  - m/sec 

RUN 

E .-j 

proj 

6 

Energy  in  Prop. 

Prop.  Mass  - kg 

RUN 

M -kg 

prop 

7 

Specific  En.-j/kg 

RUN 

avail  ■’ 

8 

Interior  losses 

chamber 

RUN 

("'’''’o'choke 

9 

Barrel  Length  - m 

RUN 

^barrel 

10 

Barrel  Dia,  - m 

RUN 

(L/D) 

gun 

11 

Groove  Hgt.  - m 

RUN 

25. 

12 

Yield  § Rad.  of 

3 

Sp.  Gravity  - kg/m 

RUN 

Zp  - m 

Eqv.  Exp. 

Card  2 

13 

Load  side  A 

Clr  2nd 

Read 

14 

Load  side  B 

2nd 

Read 

15 

Initialize 

Clr  A 

A 

16 

Overpressure 

Pressure  Ratio 

RUN 

R - m 

(P/Po) 

17 

Pulse  Length 

RUN 

Ap  - sec 

18 

Repeat  16,  17  for  futher  (P/Pq) . 
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The  following  can 

serve  as  a test  case  (M107) 

Step 

Enter 

Display 

4 

66.8  kg 

33.4 

5 

914.4  m/sec 

2.79x10^  j 

6 

25.93kg 

25.93 

7 

3.173x10^  j/kg 

5.435x10^ 

8 

3,13x103 

193.42 

9 

8,94  m 

8.94 

10 

.175 

51.1 

11 

3.66x10-3  m 

25. 

12 

1.58X10^  kg/m^ 

.7057 

13-15 

3.076x10-2 

16 

1.2 

6.044  m 

17 

5.12x10-2 

sec 

18  etc. 

Hie  following  is  a listing  of  the  program. 
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MUZZLE  BLAST  ANALYSIS  - CARD  lA 


LOG 

CODE 

KEY 

COMMENTS 

LOG 

CODE 

KEY 

000 

46 

LBL* 

040 

02 

2 

11 

A 

93 

• 

81 

HLT 

(I)  Proj  Mass 

(kg) 

01 

1 

55 

T 

04 

4 

02 

2 

95 

= 

005 

65 

X 

045 

22 

INV 

81 

HLT 

(I)  Proj  Vel 

23 

LNX 

(m/sec) 

42 

STO 

40 

X2* 

00 

0 

95 

= 

05 

5 

42 

STO 

010 

00 

0 

050 

81 

HLT 

00 

0 

E . -j 

proj 

55 

81 

HLT 

(I)  Prop  Mass 

(kg) 

81 

HLT 

42 

STO 

42 

STO 

00 

0 

00 

0 

015 

01 

1 

^prop"^8 

055 

06 

6 

65 

X 

95 

= 

81 

HLT 

(I)  Specific  Energy 

42 

STO 

(j/kg) 

00 

0 

95 

= 

07 

7 

42 

STO 

060 

81 

HLT 

020 

00 

0 

55 

■J 

02 

2 

E -j 

43 

RCL 

75 

- 

prop 

00 

0 

43 

RCL 

06 

6 

00 

0 

065 

95 

= 

025 

00 

0 

42 

STO 

95 

= 

00 

0 

42 

STO 

08 

8 

00 

0 

45 

YX 

03 

3 

E ■i-j 

070 

93 

. 

030 

81 

HLT 

avail  ■’ 

01 

1 

AP  cham 

P 

95 

= 

42 

STO 

r 

a 

20 

1/X^ 

00 

0 

65 

X 

04 

4 

AP  , 

77—  cham 

075 

01 

1 

r 

a 

05 

5 

23 

LKX 

95 

= 

035 

65 

X 

42 

STO 

93 

• 

00 

0 

09 

9 

080 

09 

9 

02 

2 

75 

- 

75 

- 

43 

RCL 

COMTvIENTS 


p—  choke 
a 

Barrel  Length-m 
Barrel  Dia-m 

Dgun"^ 

(L/D)  Gun 
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MUZZLE  BLAST  ANALYSIS  - CARD  lA  (Cont) 


LOG 

CODE 

KEY  COMMENTS 

REGISTERS 

080 

00 

0 

00 

E .-j 

07 

7 

proj  •’ 

085 

95 

= 

01 

^rop”^2 

80 

I£  Pos* 

12 

65 

B CL/D)ch>(L/D)G 

X 

02 

E -j 

prop  •’ 

090 

01 

93 

1 

03 

E • -i-j 

aval  1 

A D 

09 

9 

04 

-5—  cham 

65 

X 

r 

a 

43 

RCL 

05 

(AP/p^)  choke 

00 

0 

095 

08 

8 

06 

^barrel 

95 

= 

85 

+ 

07 

(E/D) Gun 

43 

RCL 

00 

0 

08 

h 

100 

05 

5 

23 

LNX 

09 

(L/D)  choke 

95 

= 

22 

INV 

10 

(reduction  in 

23 

LNX 

105 

20 

1/X* 

11 

^0  ■ j 

65 

X 

3 

43 

RCL 

12 

M'  - M 

00 

0 

04 

4 

13 

(1/3) 

no 

95 

= 

42 

STO 

14 

Ro  - m 

15 

Zq  - m 

16 

Z^;  - m 

^avail) 
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MUZZLE  BLAST  ANALYSIS  - CARD  IB 


LOG 

CODE 

KEY 

COMMENTS 

LOG 

CODE 

KEY 

COMMENTS 

112 

01 

1 

152 

54 

) 

00 

0 

Reduction  in  Eavail 

157 

42 

STO 

65 

X 

01 

1 

06 

6 

03 

3 

(1/3) 

95 

= 

95 

= 

117 

20 

1/X* 

42 

STO 

65 

X 

162 

01 

1 

43 

RCL 

04 

4 

Rq  - m 

00 

0 

45 

YX 

03 

3 

03 

3 

122 

95 

= 

85 

+ 

42 

STO 

167 

43 

RCL 

01 

1 

01 

1 

01 

1 

Y - j 

02 

2 

08 

8 

o 

95 

= 

127 

93 

• 

45 

yX 

02 

2 

172 

43 

RCL 

52 

EE 

01 

1 

94 

+/- 

03 

3 

03 

3 

95 

= 

132 

65 

X 

42 

STO 

43 

RCL 

177 

01 

1 

00 

0 

05 

5 

Zq  - m 

01 

1 

45 

yX 

95 

= 

03 

3 

137 

42 

STO 

95 

= 

01 

1 

182 

20 

1/X* 

02 

2 

M'  - M 

65 

X 

02 

2 

43 

RCL 

05 

5 

01 

1 

142 

65 

X 

01 

1 

81 

HLT 

(I)  Specific  Grav  - 

187 

95 

= 

kg/m^ 

45 

yx 

95 

= 

93 

• 

20 

1/X* 

03 

3 

65 

X 

00 

0 

147 

43 

RCL 

192 

07 

7 

00 

0 

07 

7 

01 

1 

65 

X 

95 

= 

01 

1 

45 

YX 

93 

• 

152 

53 

( 

197 

07 

7 

01 

1 

08 

8 

55 

52 

EE 

03 

3 

94 

+1- 
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MUZZLE  BLAST  ANALYSIS  - CARD  IB  (Cont) 


LOC 

CODE 

KEY 

02 

2 

202 

65 

X 

43 

RCL 

01 

1 

05 

5 

95 

= 

207 

42 

STO 

01 

1 

06 

6 

81 

HLT 

46 

LBL* 

212 

12 

B 

00 

0 

20 

1/X* 

81 

HLT 

COMMENTS 


t - 
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MUZZLI]  BLAST  ANALYSIS  - CARD  2A 


LOG 

CODE 

KEY 

COMMENTS 

LOC 

CODE 

KEY 

COMMENTS 

000 

46 

LBL* 

045 

00 

0 

11 

A 

00 

0 

A=QtZtq 

43 

RCL 

46 

LBL* 

01 

1 

12 

B 

06 

6 

81 

HLT 

(I)  n/Po 

005 

45 

YX 

050 

42 

STO 

03 

3 

01 

1 

75 

- 

08 

8 

43 

RCL 

65 

X 

01 

1 

06 

6 

010 

02 

2 

055 

85 

+ 

95 

=■ 

01 

1 

45 

YX 

95 

= 

43 

RCL 

55 

i 

01 

1 

53 

( 

015 

03 

3 

060 

43 

iu:i, 

95 

= 

01 

1 

42 

STO 

08 

8 

01 

1 

85 

+ 

07 

/ 

R^  - 111 

06 

6 

OJO 

01 

1 

065 

54 

) 

93 

• 

95 

= 

04 

4 

42 

STO 

20 

1/X* 

00 

0 

42 

STO 

01 

1 

1) 

025 

00 

0 

070 

75 

- 

05 

05 

5 

1/y 

01 

1 

45 

RCL 

95 

= 

01 

1 

55 

06 

6 

43 

RCL 

050 

45 

075 

00 

0 

05 

5 

01 

1 

95 

• 

95 

= 

02 

42 

STO 

05 

5 

00 

0 

055 

65 

X 

080 

02 

2 

(D-l)/D 

09 

9 

43 

RCL 

93 

• 

00 

0 

05 

5 

01 

1 

05 

5 

20 

1/X* 

040 

52 

EE 

085 

65 

X 

94 

+/- 

43 

RCL 

02 

2 

01 

1 

95 

= 

08 

8 

42 

STO 

45 

YX 

77 


CODE 

43 

00 

03 

75 

01 

95 

65 

02 

93 

05 

95 

42 

00 

04 

20 

65 

43 

00 

00 

95 

45 

93 


MUZZLE  BLAST  ANALYSIS  - CARD  2A  (Cont) 


KEY  COMMENTS 

RCL 

0 

3 

1 

X 

2 

5 

STO 

0 

4 Q 
1/X* 

X 

RCL 

0 

0 

YX 


REGISTERS 


00 

A 

01 

D 

02 

(D-l)/D 

03 

1/Y 

04 

Q 

05 

06 

07 

08 

09 

10 

11 

12 

M' 

13 

(1/5) 

14 

Ro 

15 

Zo 

16 

Zt 

17 

Rt 

18 

p/p 

a 

19 

Z 
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LOG 

112 

117 

122 

127 

132 

137 

142 

147 

152 


MUZZLE  BLAST  ANALYSIS  - CARD  2B 


CODE 

KEY 

COMMENTS 

LOC 

CODE 

KEY 

COMMENTS 

03 

3 

157 

42 

STO 

01 

1 

00 

0 

95 

= 

05 

5 

u - m/ sec 

42 

STO 

43 

RCL 

01 

1 

01 

1 

09 

9 

162 

09 

9 

45 

YX 

45 

YX 

03 

3 

53 

( 

75 

- 

93 

• 

43 

RCL 

02 

2 

01 

1 

167 

05 

5 

02 

2 

94 

+/- 

95 

= 

54 

) 

45 

YX 

65 

X 

43 

RCL 

43 

RCL 

01 

1 

172 

00 

0 

03 

3 

00 

0 

95 

= 

65 

X 

65 

X 

43 

RCL 

01 

1 

00 

0 

93 

• 

177 

02 

2 

06 

6 

65 

X 

02 

2 

05 

5 

04 

4 

55 

v 

95 

= 

Disp  R-M 

43 

RCL 

81 

HLT 

182 

01 

1 

53 

( 

08 

8 

43 

RCL 

95 

= 

01 

1 

45 

YX 

08 

8 

43 

RCL 

75 

- 

187 

01 

1 

01 

1 

03 

3 

54 

) 

95 

= 

65 

X 

55 

43 

RCL 

43 

RCL 

00 

0 

192 

00 

0 

02 

2 

05 

5 

95 

= 

95 

= 

30 

41 

GTO 

65 

X 

12 

B 

Display  At  at  B 

02  2 
08  8 
00  0 
93 

95  = 
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App-2.  Computer  Code  for  the  Theoretical  Predictions. 


The  theory  discussed  in  this  report  was  programmed  in  BASIC,  a 
listing  of  which  follows.  Execution  of  the  program  produces  output  as 
illustrated  in  Table  App  - 2.  The  initial  section  requests  the  input 
quantities  in  the  units  shown.  The  subsequent  quantities  are  calcu- 
lated and  interpreted  as  follows: 


TOTAL  ENERGY  IN  PROPELLANT 

Total  energy  available 

KINETIC  ENERGY  IN  PROJECTILE 

Projectile  KE  at  the 
nominal  muzzle  velocity 

MAXIMUM  AVAILABLE  ENERGY 

"t  ■ '"ke’ 

OVERPRESSURE  RATIO  (CHOKE) 

Results  of  Bernoulli 
expansion 

OVERPRESSURE  RATIO  (MUZZLE) 

Includes  energy  impedance 
from  rough  wall 

YIELD 

Yield  of  spherical  explosion 

SPHERICAL  RADIUS 

Radius  of  sph.  explosion 

TRANSITION  RADIUS  (R^) 

Radius  separating  strong  and 
weak  shock  regimes 

TRANSITION  VELOCITY  (U^) 

Material  speed  at  R^ 

The  remaining  tabulation  gives  the  range,  predicted  overpressure  ratio, 
pulse  length  and  comparisons  between  the  predicted  results  and  the  noise 
level  restrictions  specified  by  the  safety  regulations^ ^ . The  last  two 
columns  are  interpreted  as  the  number  of  shots  per  day  which  can  be 
tolerated  with  single  or  double  ear  protection. 


O -7 

"Military  Standard  - Noise  Limits  for  Army  Material,"  MIL  - STD  - 
14?4A(MI,  2 Mar.,  1975. 


80 


Table  App-2.  Muzzle  Blast  Overpressure 
and  Pulse  Length 


PROJECTILE  MASS  (KG) 

VELOCITY  (M/SEC) 


? 66.8 

? 914.4 


PROPELLANT  MASS  (KG) 

SPECIFIC  GRAVITY  (KG/M''3) 
SPECIFIC  ENERGY  (J/KG) 


? 25.93 

? 1.58E3 

? 3.173E6 


CHAMBER  PRESSURE  RATIO  (P/PO) 


3.13E3 


BARREL  LENGTH  (M) 

DIAMETER  (M) 
GROOVE  HEIGHT  (M) 


? 8.94 

? .175 

? 3.66E-3 


ENERGY  CALCULATIONS 


TOTAL  ENERGY  IN  PROPELLANT 
KINETIC  ENERGY  IN  PROJECTILE 
MAXIMUM  AVAILABLE  ENERGY 


8.22759E+7  JOULES 
2.79263E+7  JOULES 
5.43496E+7  JOULES 


INTERIOR  LOSSES 

OVERPRESSURE  RATIO  AT  3.86443  M IS  195.476  (CHOKE) 

OVERPRESSURE  AT  8.94  M IS  62.7723  (MUZZLE) 

EQUIVALENT  EXPLOSION  PARAMETERS 


YIELD  181722  JOULES 

SPHERICAL  RADIUS  8.67565E-2  M 


MUZZLE  BLAST  CALCULATIONS 


TRANSITION  RADIUS  (RT) 
TRANSITION  VELOCITY  (UT) 


SPHERICAL 

CYLINDRICAL 


.543517  M 
.882521  M 
172.67  M/SEC 
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Table  App-2.  Muzzle  Blast  Overpressure 
and  Pulse  Length  (Cent) 


RADIUS 

OVERPRESSURE 

PULSE  LENGTH 

# ROUNDS 

W/ PROTECTION 

CM) 

(P-P0)/P0 

T/ (RT/UT) 

SINGLE 

DOUBLE 

.792133 

.410223 

NONE 

NONE 

4 

.328743 

.736885 

5 

100 

6 

.201607 

1.00585 

5 

100 

8 

.143676 

1.24237 

100 

1000 

10 

.112007 

1.44735 

100 

1000 

12 

9.17494E-2 

1.63398 

100 

1000 

14 

7.81999E-2 

1.7987 

100 

1000 

16 

6.92267E-2 

1.93242 

1000 

>1000 

18 

6.21047E-2 

2.05978 

1000 

>1000 

20 

5.73473E-2 

2.15503 

1000 

>1000 
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MUZBLA 


100  FOR  N=1  TO  5 STEP  1 
110  PRINT 

120  NEXT  N 

140  PRINT  "MUZZLE  BLAST  OVERPRESSURE  AND  PULSE  LENGTH" 

150  PRINT  " " 

160  PRINT 

170  PRINT 

200  REM  ....  INPUT  SECTION  .... 

201  REM  . . VARIABLE  DEFINITIONS  ARE  AS  FOLLOWS: 

202  REM  . . PROJECTILE:  VO  - VELOCITY  (M/SEC),  WO  - MASS  (KG) 

203  REM  . . PROPELLANT:  W1  - MASS  (KG),  XI  - SP.  GRAVITY  (KG/M'S) 

204  REM  . . X2  - SPECIFIC  ENERGY  (J/KG) 

205  REM  . . PRESSURE:  P - CHAMBER  PRESSURE  RATIO  (P/PA) 

206  REM  . . BARREL:  LO  - Length  (M) , DO  - Diameter  (M) , 

207  REM  . . LI  - GROOVE  HEIGHT  (RIFLING)  (M) 

210  PRINT  "PROJECTILE",  "MASS  (KG)",, 

220  INPUT  WO 

225  PRINT 

230  PRINT  , "VELOCITY  (M/SEC)", 

240  INPUT  VO 

245  PRINT 

246  PRINT 

250  PRINT  "PROPELLANT", "MASS  (KG)",, 

260  INPUT  W1 

265  PRINT 

270  PRINT  , "SPECIFIC  GRAVITY  (KG/M''3)", 

280  INPUT  XI 

285  PRINT 

290  PRINT. , "SPECIFIC  ENERGY  (J/KG)", 

300  INPUT  X2 

305  PRINT 

306  PRINT 

310  PRINT  "CHAMBER  PRESSURE  RATIO  (P/PA)",, 

320  INPUT  PO 

325  PRINT 

326  PRINT 

330  PRINT  "BARREL", "LENGTH  (M)",, 

340  INPUT  LO 

345  PRINT 

350  PRINT  , "DIAMETER  (M)",, 

360  INPUT  DO 

365  PRINT 

370  PRINT  , "GROOVE  HEIGHT  (M)",. 

380  INPUT  LI 

390  PRINT 

395  PRINT 

400  REM  . . END  INPUT  DATA 
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1000 

1010 

1011 

1020 

1030 
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REM  . . 

REM  . . BEGIN  ENERGY  CALCULATIONS 

REM  . . E(PROP)  = E2,  E(PROJ)  = El,  E (AVAIL)  = EO 

LET  E2=X2*W1 

LET  E1  = .5*W0*V0''2 

LET  E0=E2-E1 

PRINT 

PRINT  ENERGY  CALCULATIONS" 

PRINT 

PRINT  "TOTAL  ENERGY  IN  PROPELLANT" , ,E2 ;" JOULES" 

PRINT  "KINETIC  ENERGY  IN  PROJECTILE", , El ; "JOULES" 

PRINT  " MAXIMUM  AVAILABLE  ENERGY", ,E0;"  JOULES" 

LET  P0=P0-1 
PRINT 

PRINT  ,"  INTERIOR  LOSSES" 

PRINT 

REM  . . ROUGHNESS  FACTOR  = H,  OVERPRESSURE  RATIO  AT  CHOKE  = X 

REM  . . O'PRESSURE  RATIO  AT  MUZZLE  AND  INITIAL  YIELD  = YO 

REM  . . CHOKE  LENGTH  = L2 

LET  . . LET  H=L1/D0 

LET  L2=0 

LET  X=P0 

REM  . . IF  CHOKE  L/D  > GUN  L/D  ASSUME  NO  CHOKE  LOSSES 
IF  H=  0 GOTO  2180 
IF  15/H''.1>(L0/D0)  GOTO  2180 
lET  L2=15/H^.l 

REM  . . IF  P0>100  USE  FITTED  CURVE 
IF  P0>100  GOTO  2160 
LET  X=P0 

LET  C= (1+X) * (l+4*X"2/5/ (1+X) / (9+X) ) "5-1 
IF  ABS  (C-P0)<.001  GOTO  2170 
LET  X=X-(X-1)*(C-P0)/C-1.29) 

GOTO  2110 

LET  X=  EXP  (.9211*  LOG  (P0)-2.138) 

PRINT  "OVERPRESSURE  RATIO  AT";L2*D0;"  M IS",X;"  (CHOKE)" 

LET  Y0=  EXP  ((2*H*(L2-L0/D0)+1.068*  LOG  (X))/1.068) 

PRINT  "OVERPRESSURE  AT";L0;"  M IS",,Y0;" 

LET  F=PO/YO 

REM  . . REDUCTION  IN  E AVAILABLE 
LET  YO=EO/(6*F) 

PRINT 

PRINT  ,"  EQUIVALENT  EXPLOSION  PARAMETERS" 

PRINT 

PRINT  "YIELD",,, YO'"  JOULES" 

LET  R0= (Wl/ (4*3. 14159/3*6*X1)) " (1/3) 

PRINT  "SPHERICAL  RADIUS" ,, RO; " M" 

REM  . . END  INTERIOR  LOSS  CALCULATIONS 
REM  . . 
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3001  REM  . . BEGIN  OVERPRESSURE  VS,  DISTANCE  CALCULATIONS 

3002  REM  . . 

3010  REM  , , ATMOSPHERIC  CALCULATIONS  AND  DEFINITIONS 
3020  REM  . . PRESSURE  (J/M'^S)  = P5,  DENSITY  (KG/M''3)  = D5 
3030  LET  P5=101325 
3040  LET  D5=1.29 

3050  REM  . . EXPONENTS  Q1§Q2,  M IN  MEZ,  G=GAMMA,  Q9=Q(TRAN) , Z9=Z (TRAN) 
3060  LET  Ql=3.25 
3070  LET  Q2=3.25 
3080  LET  G=1.4 

3090  LET  M=.25*W1/ (4/3*3.14159*6*05) 

3100  LET  Z0=(R0"3+M)^(l/3) 

3110  LET  Q9=.0955 

3120  REM  . . SOLVE  FOR  TRANSITION  RADIUS 
3130  REM  . . GET  CORRECT  DIMENSIONLESS  YO 
3140  LET  Y0=.00001*Y0 
3150  LET  A=Y0/(4*3.14159*Q9*Z0''3) 

3160  IF  Q1=Q2  GOTO  3240 
3170  LET  B=(Q2-Ql)/(Q2-3)/(Ql-3) 

3180  LET  V=4 

3190  LET  T=(V"Ql)/(Ql-3)-B*V''3 

3200  IF  ABS  (T/A-l)<.001  GOTO  3250 

3210  LET  Vl=(A-T)/(Ql*(V"(Ql-l))/(Ql-3)-3*V''2*B) 

3220  LET  V=V+V1 
3230  GOTO  3190 
3240  LET  V=(A/4)"(1/Q1) 

3250  LET  Z9=V*Z0 

3251  LET  R9=(Z9"3-M)"(l/3) 

3252  LET  R8=. 8735* (1 . 5* (Ql-2)/ (Ql-3) ) ^ (1/Ql ) 

3253  REM  . . NOTE:  R8  IS  CYL.  ADJUSTMENT  ASSUMING  Q1=Q2 

3254  REM  . . THERE  IS  NOT  INCLUDED  CALC.  FOR  UNEQUAL  Q'S 
3260  PRINT 

3270  PRINT,”  MUZZLE  BLAST  CALCULATIONS” 

3280  PRINT 

3290  PRINT  "TRANSITION  RADIUS  (RT) ”, "SPHERICAL”,  R9;”  M” 

3291  PRINT  ,, "CYLINDRICAL”, R8*R9;”  M” 

3300  LET  U9=172.67 

3310  PRINT  "TRANSITION  VELOCITY  (UT)”,,U9;”  M/SEC” 

4000  REM  . . BEGIN  QZQ  CALCULATIONS 

4001  REM  . . A = Q(T)*Z(T)"Q1,  B = Q(T) *Z (T) "Q2 , Q = WASTE  HEAT 

4002  REM  . . PRINT  HEADER 

4003  PRINT 

4004  PRINT  "RADIUS”,  "OVERPRESSURE” , "PULSE  LENGTH” , ”#ROUNDS  W/  PROTECTION" 

PROTECTION 

4005  PRINT  " (M)”,”  (P-P0)/P0”,”  T/ (RT/UT)”, "SINGLE”, "DOUBLE” 

4006  PRINT  " "," "," " ” " 

4009  LET  P=2 

4010  LET  R=2/R8 
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IF  R>R0  GOTO  4050 
LET  R=R+2/R8 
GOTO  4020 
LET  Z=(R^3+M)"(l/3) 

LET  Q5=Q1 
IF  Z<Z9  GOTO  4090 
LET  Q5=Q2 

LET  Q=Q9*(Z9/Z)"Q5 

IF  Q<2.718  GOTO  4156 

REM  . . HIGH  PRESSURE  APPROXIMATION 

LET  X=.5*(23-  SQR  (441-64* .4343*  LOG  (Q))) 

LET  P=10^X+1 
GOTO  4220 

REM  . . MEDIUM  PRESSURE  Q - NFWTON/RAPHSON  ITERATION 
LET  A1=(G+1)/(G-1)*P+1 
LET  A2=P+(G+1)/(G-1) 

LET  F=(A2*P^(1/G)-A1)-A1*(G-1)*Q 
LET  F1=1/G*P^(1/G-1)*A2+P^(1/G) 

LET  F1=F1-(G+1)/(G-1)-(G+1)*Q 
LET  N=F/F1 

IF  ABS  (N)<.001  GOTO  4220 
LET  P=P-N 
GOTO  4156 

REM  . . PRESSURE  RATIO  OBTAINED,  NEXT  DELTA  T 
REM  . . SECTION  FOR  DELTA  T 

REM  . . Y = PROMPT  ENERGY,  Y = ADJUSTED  ENERGY 

LET  A=Q9*Z9''Q1 

IF  Q1=Q2  GOTO  5130 

IF  Z<Z9  GOTO  5060 

LET  A=0 

LET  Z5=Z9 

IF  Z<Z9  GOTO  5090 

LET  Z5=Z 

LET  B2=A/(3-Ql)*(Z5"(3-Ql)-Z"(3-Ql)) 

LET  B=Q9*Z9"Q2 

LET  Y5=4*3.14159*(B2-B/(3-Q2)*Z5"(3-Q2)) 

GOTO  5140 

LET  Y5=4*3. 14159*A*Z" (3-Ql)/ (Ql-3) 

LET  D=((G+1)/(G-1)*P+1)/(P+(G+1)/(G-1)) 

LET  Y=Y5/ (4/3*3. 14159*P)/(1+G)*(D-1)/D 
LET  U=  SQR  (P5/D5)*  SQR  ((P-1)* (D-l)/D) 

LET  T=Y"(l/3)/U 

REM  . . OUTPUT  MUZZLE  BLAST  VALUES 
PRINT  R*R8, (P-1) ,T/(R8*R9/U9) , 

REM  . . THIS  SECTION  REFLECTS  THE  SG'S  SAFETY  REQUIREMENTS 
REM  . . O' PRESSURE  AND  PULSE  LENGTH  COMPARED  TO  STANDARD 
LET  H9=-3.095*  LOG  (T*1000) + 189 . 9 
LET  H8=. 4353*20*  LOG  ( (P-1) *14 . 7*3 .4475E+8) 
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6040  IF  H8>H9  GOTO  6170 
6060  IF  H8>(H9-6.5)  GOTO  6150 
6070  IF  H8>(H9-11.5)  GOTO  6130 
6080  IF  H8>140  GOTO  6110 
6090  PRINT  ">1000", ">1000" 

6100  GOTO  7000 

6110  PRINT  " 1000", ">1000" 

6120  GOTO  7000 

6130  PRINT  " 100","  1000" 

6140  GOTO  7000 

6150  PRINT  " 5","  100" 

6160  GOTO  7000 

6170  PRINT  " NONE","  NONE" 

7000  LET  R=R+2/R8 

7010  IF  R*R8>21  GOTO  9000 

7020  GOTO  4050 

9000  PRINT 

9010  PRINT 

9020  PRINT  "OVERPRESSURE  (DB)  AT  ";R*R8;"  M IS";H8 
9030  END 
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A 

B 

D 

E 

F 

H 

I 

K 

L 

M,  M' 

N 

P 

Pr 

APr 

Q 

R 

S 

T 

U 

V 
W 
X 

Y 
Z 


Cross-sectional  area  of  barrel 

Ratio  of  mass  prompt  energy  to  that  of  air; 
MEZ(II-B.2) 

Barrel  hydraulic  diameter 
Energy 

Gibbs  free  energy 

Ratio  of  roughness  height  to  tube  diameter,  h/D 
Function  in  Equation  (17) 

Kinetic  energy  per  unit  mass,  1/2  u^ 

Length  (fixed)  along  barrel  axis 
Mass,  Mass  function  (MEZ) 

Number  of  Particles  (kinetic  theory) 

Absolute  pressure 
Pressure  ratio,  P/Pa 
Overpressure  ratio,  (P-Pa)/Pa 
Waste  heat 

Distance  from  point  of  explosion  to  shock  wave 

Perimeter  of  barrel 

Temperature 

Interaction  potential 

Speed,  volume 

Pressure  volume  work,  W = /Pdv 
Dimensionless  distance  (Z  /Z  ) 

Prompt  Energy 

Mass  corrected  R,  MEZ 
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a 

e 

h 

k 

q 

r 

t 

u 

V 

X 

z 


Speed  of  sound 
Specific  energy,  (E/mass) 

Height  of  barrel  roughness 
Boltzmann's  constant 
Exponent  in  QZQ  hypothesis 
Radial  coordinate 
Time 

Material  velocity  behind  shockwave 
Specific  volume 

Dimensionless  distance,  (L/D) , in  gun 
Coordinate  along  barrel  axis 


Subscripts 
a Ambient 

c Choke  joint 

f Final 

I Initial,  internal 

K.E.  Kinetic  energy 
p Peak  value 

T Total 

o At  charge 

t Transition 


Greek 

a Proportionality  constant 

6 2y/(y-1) 

Y Ratio  of  specific  heats 

e Energy  per  unit  volume 

5 n Constants  in  intermediate 

calculations 

y (y+i)/(y-i) 

p Density 

6 Angle  between  unit  normal 
along  barrel  axis  and  posi- 
tion vector  from  muzzle  to 
pressure  probe. 
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